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y—( ■ Abstract. We review the status of the two-dimensional Abelian sandpile model as a strong 

candidate to provide a lattice realization of logarithmic conformal invariance with central 
charge c = —2. Evidence supporting this view is collected from various aspects of the 
^ ■ model. These include the study of some conformally invariant boundary conditions, and the 

| corresponding boundary condition changing fields, the calculation of correlations of certain 

bulk and boundary observables (the height variables) as well as a proper account of the 
necessary dissipation, which allows for a physical understanding of some of the strange but 
generic features of logarithmic theories. 
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1. Introduction 

The Abelian sandpile model was originally conceived as a prototypical example of a so-called 
self-organized critical system. Namely, a dynamical system which automatically adjusts its 
own behaviour to converge and maintain itself in a critical state. According to its inventors, 
Bak, Tang and Wiesenfeld [1], the sandpile model is just a specific incarnation, in a concrete 
situation, of a generic mechanism present in Nature, which would explain why so many power 
laws are actually observed (earthquakes, avalanches, solar flares, ...). 

Although criticality was built in the model from the start, the focus was rather on 
dynamical aspects and the extent to which the dynamics defining the sandpile model is 
universal. It was also soon realized that the asymptotic regime is described by a critical 
measure, and therefore lends itself to the usual analyses, like scaling and conformal invariance, 
especially in two dimensions. The first hint of the central charge (c = —2) as well as the 
non-local features of the model led to the suspicion that the conformal description could very 
well be logarithmic. 

If the abstract understanding of the general features of logarithmic theories has greatly 
improved, as this Special Issue clearly shows, the situation is not yet satisfactory in most 
cases. In this regard, lattice models realizing logarithmic conformal invariance may prove 
extremely useful. Infinitely many examples are known by now, some being better understood 
than others, but for most of them, it turns to be extremely difficult to completely pinpoint 
the specific conformal theory at work, a rather unusual situation if we compare with the 
unitary conformal theories. 

This review will focus exclusively on one of these lattice models, the two-dimensional 
Abelian sandpile model (ASM). In most lattice realizations, the conformal spectrum is more 
easily accessible while correlators are harder, and in many cases, unknown. The ASM is 
one of the very few models where the situation is the inverse: correlation functions of 
some observables can be computed exactly, in the bulk, on boundaries, or both, whereas 
the complete conformal spectrum is not easy to determine because a transfer matrix 
formulation is not always available. This peculiarity makes the ASM rather unique among 
the lattice realizations of logarithmic conformal invariance, and helps to develop a physical 
understanding of this class of conformal theories. 

We will start by giving a general overview of the ASM, as a discrete dynamical system 
in 2 + 1 dimensions, and its relationship to other models (spanning trees, dimers, loop-erased 
random walks). In the asymptotic, stationary regime, the statistical properties of the model 
are controlled by a probability measure P on the space of configurations of the sandpile. For 
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a finite system, this measure is simple but non-local with respect to the natural degrees of 
freedom, namely the heights of the pile at each site, which are strongly correlated throughout 
the lattice. Such intrinsically non-local features are a landmark of LCFTs, and qualifies the 
ASM as a good candidate for an LCFT description. Indeed most explicit computations to 
date aim at providing evidence that the scaling limit of the measure P is the field-theoretical 
measure of an LCFT, with central charge c = — 2. 

The rest of the article will be devoted to review the exact results which are the most 
convincing to support the previous assertion. Such calculations include the calculation of 
certain bulk correlators, the discussion of a few boundary conditions and the corresponding 
boundary condition changing fields, and some boundary and bulk correlators, some of them 
forming Jordan cells. In the ASM, an absolutely crucial role is played by dissipation, which 
is essential to make the dynamics well-defined. In the conformal picture, the insertion of 
dissipation is represented by a dimension field, logarithmic partner of the identity, which 
allows for a transparent understanding of some of the strangest features of LCFT. 

We should mention that the ASM is not the only lattice model believed to be described 
by a logarithmic conformal theory with c = —2. At least two other models are known, namely 
the dense polymer model which, in its loop formulation, is the first of the infinite series of 
so-called logarithmic minimal models [2], and the dimer model [3]. In their scaling limit, the 
three models are believed to be different. Although very close and even equivalent in certain 
instances, the ASM and the dimer model are not when the lattice has one or more periodic 
direction (see below in Section 2.4). Some differences between the dense polymer model and 
the dimer model have been recently emphasized in [4]. Finally, while the dense polymer is 
likely to be described by the symplectic free fermion theory, it is not the case of the ASM, 
as an argument recalled in Section 5.1 shows. Exactly which conformal theories describe the 
dimer model and the ASM is a widely open question. 

2. The Abelian sandpile model 

We review in this Section the definition of the model and its basic features, omitting most of 
the time the detailed proofs. These and more details about the model can be found f.i. in 
the reviews [5, 6]. We also point out its relation to other mathematical problems. 

2.1. Definition 

The sandpile model we are going to consider has been defined by Bak, Tang and Wiesenfeld 
in 1987 [1], but the first systematic mathematical analysis has been carried out by Dhar [7], 
who showed in particular its Abelian structure. The BTW model was subsequently called 
the Abelian sandpile model (ASM). The model was originally defined on a finite portion of 
Z d , but can as easily be defined on any finite graph C, chosen unoriented to avoid technical 
subtleties [8]. As we will be primarily concerned with two-dimensional lattices, we will mostly 
think of C as being a finite portion C of Z 2 , typically a rectangular LxM grid, with iV = LM 
sites. The formulation and the results which follow are however general. 

It is convenient to consider an extended graph £*, which is the graph £ itself, 
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supplemented by one extra vertex called the root and edges connecting the root to sites 
of C. For each vertex i, we accordingly define z i as the coordination number of % in C and z* 
as its coordination number in £*, so that z\ — z i > is the multiplicity of the edge connecting 
i to the root. We also introduce the symmetric toppling matrix A on £ as follows, 

{z* for i = j, 
— 1 if % and j are neighbours, (2.1) 
otherwise. 

The ASM is a stochastic dynamical open system, discrete in space and in time. The 
microscopic degrees of freedom are random variables hi, called heights, attached to the sites: 
hi counts the number of sand grains at site % and takes integer values, conventionally chosen 
to be larger or equal to 1. A configuration is a set of values {hi}i € c; it is called stable if all 
heights satisfy 1 < hi < z*. We denote by S the set of all stable configurations, of cardinal 

The dynamics is defined on S as follows. Given a configuration C t at time t, the dynamics 
produces the configuration C t+1 in two steps. 

(a) Seeding. A grain is dropped on a random site j according to a distribution {pj} on 
C, that is, hj — > hj + 1, and hi — > hi for all i ^ j. C t +\ is this new configuration if it 
is stable; if it is not stable (new hj = z* + 1), we go to step (b). 

(b) Relaxation. Every unstable site -with a height larger or equal to z* + 1- topples: it 
loses z* grains of sand and gives one to each of its Zj nearest neighbours in C. In other 
words when a site j topples, the whole configuration is updated as hi — > hi — Ajj. 
When all unstable sites have toppled, the resulting configuration is stable and defines 
Ct+i- 

We note that after a toppling at j, the sand is redistributed but remains within C if 
j is not connected to the root (z* = z-); if j is connected to the root, z* — Zj grains are 
transfered to the root and never come back in C. In the former case, we say that the site j 
is conservative, while in the latter case, it is called dissipative. 

The above dynamics is well-defined. The seeding of a grain at j can potentially trigger 
a large avalanche of topplings. If unstable after step (a), the site j topples and transfers 
sand to its neighbours, which can themselves become unstable and topple, making their own 
neighbours unstable, and so on. However the relaxation process always terminates provided 
the set of dissipative sites is not empty and every site is path-connected to at least one 
dissipative site. Indeed if the avalanche were to propagate for ever, an infinite number of 
grains would be evacuated to the root, which is impossible since any stable configuration 
holds a finite quantity of sand. Thus the root acts like a infinitely deep sink to which sand 
is irreversibly lost. The root is therefore crucial to make the dynamics well-definedf . 

In the cases that are most usually considered, namely a rectangular grid or a cylindrical 
grid, the dissipative sites are chosen among the boundary sites. A boundary site connected 
to the root is called open, and closed if it is not connected to the root. Similarly a whole 

f This would not be the case if the model was defined from the start on an infinite lattice. The set of allowed 
configurations must then be appropriately restricted [6]. 
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boundary consisting of open resp. closed sites is called an open resp. closed boundary. The 
bulk sites are usually not connected to the root, and so have z i = z* = 4 on the square 
lattice. The toppling matrix becomes a Laplacian-like operator with appropriate boundary 
conditions. 

The relaxation process is also independent of the order in which the topplings take place. 
This is a consequence of the fact that whether or not a toppling occurs at site j depends only 
on hj, and that the effects of two topplings at j and k commute, 

hi^hi- Aij - A ijfc = h t - A ijfc - Aij. (2.2) 

If we denote by a,j the operator that drops a sand grain at j and executes the relaxation 
process, we can rephrase the dynamics by saying that given C t , it returns Ct+i = ctjC t with 
probability pj. The operators a,-, for j e C, all map the set of stable configurations into itself, 
and from the previous observations, are mutually commuting, 

[a l ,a J }=0, Vi,je£. (2.3) 

The qualification of this sandpile model as Abelian stems from this property. 

The dynamics as defined above allows to follow the stochastic time evolution of any 
given initial configuration. As we are interested in the long-time statistical properties of 
the configurations Ct, it is more convenient to compute the time evolution of probability 
distributions P t (C) over the set of stable configurations. From an initial distribution Po(C), 
one easily obtains the following master equation, 

P t (C) = Pi E S ( C - a > C ') ^-i(C) = E W W P *-i( C ')- (2-4) 
jec c c 

The matrix Wc,c is the transition matrix for a Markov chain with finite state space S, 

equivalently of a random walk on S. 

The long-time behaviour of the sandpile is controlled by the limit measure ¥* c = 
lim^oo P t = lim^oo W t F . The general theory of Markov chains applied to this particular 
case yields the following results. 

The set of stable configurations is divided into two disjoint subsets: the recurrent 
configurations are in the repeated image of the dynamics and keep reoccuring infinitely 
often; those which can only occur a finite number of times are called transient. For all t large 
enough, P t (C) = if C is transient and P*(C) > if C is recurrent. Therefore if we start from 
any stable configuration and apply the dynamics, it takes a finite time (of order iV 2 ) to reach 
a recurrent configuration, after which all subsequent configurations are recurrent (from the 
definitions, a recurrent configuration can never be mapped to a transient one). 

Under the hypothesis that pj > for all j G C (the uniform distribution pj — jj is 
usually used), one can show that the limit P^ is unique, i.e. does not depend on the initial 
distribution, and moreover is uniform on the set 1Z of recurrent configurations, 

{— — if C is recurrent, 
\n\ (2.5) 
if C is transient. 

The stationary distribution F* c for a given graph looks very simple, since it is uniform, 
but at the same time is extremely complex because its support, the set 1Z, has a complicated 
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structure, as we will see below. In fact, a configuration must satisfy non-local contraints to 
be recurrent. 

In the infinite volume limit, the sequence of measures P^ converges to a measure on an 
infinite lattice, which we simply denote by P. The main purpose of this review is to give 
support to the statement that 

The scaling limit of the measure P is the field-theoretic measure of a two-dimensional 
logarithmic conf or mal field theory with central charge c = — 2. 

Evidence will be collected from results obtained during the last ten years. 
2.2. Recurrent configurations and spanning trees 

The recurrent configurations manifestly play a central role. The few definitions recalled above 
allow us to gain some insight into what they really are. We assume all pi > 0. 

A first observation is that the full (or maximal) configuration C* with hi = z* is 
certainly recurrent. Indeed it can be obtained from any stable configuration by applying 
a suitable sequence of a; (by filling the sites which are strictly below their threshold). 
Then any configuration of the form Y[j a /C* obtained from C* by applying any string of 
a,j is also recurrent. Conversely, any recurrent configuration must be of this form, because 
once C* is produced by the dynamics, they are the only ones to ever appear. As any 
recurrent configuration can be obtained from C* and vice- versa, it follows that any recurrent 
configuration can be obtained from any other by applying a suitable string of a«; they form 
under the dynamics a single irreducible component. Another simple corollary is that if 
C = {hi} is recurrent, then any stable configuration C with heights h\ > hi is also recurrent. 

The previous observations imply that the operators are invertible on 1Z. For every i, 
an operator fej, monomial in the a's, should exist such that biOiC* = a^iC* = C*. It must 
also be unique because the only recurrent configuration d such that a^Cj = C* has heights 
equal to h k — z\ — S^i- Thus the operators a, restricted to the recurrent configurations have 
a unique inverse and generate a finite Abelian group. 

The question nonetheless remains: are all stable configurations recurrent ? The negative 
answer stems from the toppling rules. 

To take a simple example, consider a stable configuration in which two neighbouring 
sites have a height equal to 1. Assume also that the dynamics has been run for long enough 
so that these two sites have both toppled. The site that has toppled last has given one 
grain to its neighbour which therefore must have a height at least equal to 2. Thus any 
stable configuration with two neighbouring l's cannot be recurrent. The same argument 
applies to three neighbouring sites with height values 121 (the 2 in central position) or to 
four neighbouring sites forming a square, all with a height 2. 

More generally, a sub configuration on a subset of sites F is forbidden in a recurrent 
configuration if each site of F has a height which is smaller or equal to the number of its 
neighbours in F. The three examples given above are the smallest forbidden clusters; however 
the size of a forbidden cluster is only bounded by the size of C. For instance the configuration 
hi = Zi on C is forbidden for F = C and contains no smaller forbidden sub configuration. 
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Forbidden sub configurations form the key concept because one can prove that a 
configuration is recurrent if and only if it contains no forbidden subconfiguration [9]. Because 
one may have to scan a large portion of the lattice to see if it contains a forbidden 
subconfiguration, the recurrence property is non-local. This is what makes the sandpile 
model non-trivial, interesting but also notoriously difficult to handle analytically t.. 

The burning algorithm [7] provides a simple and practical way to test whether a 
configuration is recurrent and at the same time establishes [9] a bijective correspondence 
with rooted spanning trees on £* (hence an equivalence of the sandpile model and the q — > 
limit of the q- Potts model). It goes like follows. 

To start, let /C = £ be the set of unburnt sites at time 0. The sites of JC whose height 
is strictly larger than their number of neighbours in /Co are called burnable at time 1 (these 
sites must be dissipative). We burn them, and doing so we define a new set of unburnt sites 
/Ci Q K>o at time 1. Next we burn the sites of /Q which are burnable at time 2 (those with 
a height larger than their number of neighbours on /Ci) and obtain yet a new set /C2 Q /Ci 
of unburnt sites at time 2. This procedure is carried on until, for some T, no site of /Ct is 
burnable, implying /C^+i = fcr- If = {0}) the configuration is recurrent; otherwise, the 
subconfiguration on Kt is forbidden and the configuration is transient. 

In the sequence of burnings as described above, it is not difficult to see that a site j of /Q 
which is burnable at time t + 1 has a least one neighbour that was burnt at time t (otherwise 
j would have been burnable at an earlier time). If j has exactly one such neighbour k, we 
say that j catches fire from k or that the fire propagates along the bond {k,j). If there are 
more than one, we use an ordering prescription to decide which neighbour k sets j afire. 
The prescription is arbitrary -and may change from site to site- but it must be fixed, i.e. 
identical for all configurations. We also conventionally decide that the sites of /Co which are 
burnable at time 1 catch fire from the root. If the configuration is recurrent, the fire path is a 
collection of edges which contains no loop and such that every site is the endpoint of at least 
one edge. In other words, the fire path forms a rooted spanning tree, growing from the root 
towards the interior of C Conversely, and given a fixed ordering prescription, every spanning 
tree is associated with a unique recurrent configuration. Therefore the measure induced on 
the spanning trees is also uniform. 

Note that the rooted spanning tree one obtains from this algorithm is connected on 
£*. If one removes the root, and therefore focus on C, the tree in general will appear to be 
disconnected, i.e. will form a spanning forest. 

We can observe that the rooted spanning tree carries a natural orientation, defined by 
the fire propagation. It is however conventional to orient the rooted spanning tree in the 
opposite direction, so that the oriented edges are like arrows that flow towards the root. A 
recurrent configuration can then be viewed as a configuration of arrows: each site has one 
outgoing arrow pointing to one of its neighbours in C* (it can be the root for dissipative 
sites) with the global constraint that the arrows do not form loops. 

The total number of arrow configurations is clearly equal to Yli^c z t = l"^l> but ^ ne 

I Notable exceptions are Bethe lattices [10] or ld-like lattices [11, 12], but admittedly they show less 
interesting behaviour. 
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number of those without loops is smaller. By the matrix-tree (Kirchhoff) theorem, this 
number, also equal to the number of recurrent configurations, is equal to 

\K\ = det A. (2.6) 

We note that if C contains no dissipative site, that is z* = z i for all i, det A vanishes and no 
configuration is recurrent (there are no spanning trees). 

For C a thick domain in Z 2 containing iV sites in which the interior sites are conservative, 
one can show that the number of recurrent configurations on £ grows exponentially like 
\1Z\ ~ e^^ ~ (3.21) N with G = 0.915965 the Catalan contant, to be compared with the 
number ~ 4^ of all stable configurations. 

Thus the ASM configurations can be described in terms of height variables or spanning 
trees. The heights may look more natural, they are local but globally constrained. The global 
constraints induced by the recurrence property are encoded into the global structure of the 
spanning trees. In actual computations, the latter often prove more convenient. 



2.3. The sandpile group 

The ASM offers a large number of interesting and challenging mathematical problems [13]. 
A fine example (see Section 4.2 for others) is the sandpile group, also called the critical group 
of a graph. 

We have shown that the operators generate a finite Abelian group G (it depends on C 
but we omit the explicit dependence). The order of G is equal to \R\, because all recurrent 
configurations can be obtained from any fixed one by repeated applications of the aj. The 
recurrent configurations themselves form an Abelian group isomorphic to G, with a group 
law given by the sitewise addition of heights (followed by relaxation) 

C®C = \{a h SC' = \{afc. (2.7) 

In particular, there must be a unique recurrent configuration which is the identity in G. The 
resolution of G into cyclic factors and the geometric structure of the identity configuration 
are challenging questions. 

The addition of A^ = z* grains at site i of a recurrent configuration makes it topple. 
After the toppling of i, the height hi returns to its original value while the neighbouring 
heights have been incremented by 1. Therefore adding A iti grains at i is the same as adding 
one grain to each neighbour of i. This yields the following relations on the generators of G, 

J] af' j = 1, Vi G C. (2.8) 

jec 

Being finite Abelian, the irreducible representations of G are all one-dimensional, given 
by dj = e 2mVj . The relations (2.8) imply that £V Aijipj = rrii for an integer vector rh G Z N , 
and so (pj = Ylkec (^)jfe mfc are rational. Since two vectors m and m + An yield identical 
representations (same phases ipj), the dual group of G and thus G itself is given by 

G = Z N /AZ N . (2.9) 

In particular we recover the order \G\ — det A. 
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In terms of recurrent configurations, this last equation has an interesting consequence. 
Let us consider the set of all configurations with heights hi G Z and the equivalence C ~ C if 
and only if C and C are related by topplings, i.e. C — C + An for some integer n. Then the 
quotient Z* N /AZj N is the number of equivalence classes, and (2.9) shows that there is exactly 
one recurrent configuration in each class. In this respect, toppling invariant functions on 1Z 
are particularly useful as coordinates on G. 

The group structure of G is identical to Z^/A Z N as a finite Z-module. On general 
grounds [14], it is a direct product of finite cyclic groups 

G = Z dl x Z d2 x . . . x Z dN , d i+l \di. (2.10) 

The integers di, called the elementary divisors of A, are related to the Smith normal form§ of 
A (generically, most of them are equal to 1). For a specific, small graph, the computation of 
the elementary divisors is standard, but to compute them for large graphs or infinite classes 
of graphs is highly non-trivial. The generators of the cyclic factors in (2.10) as well as a 
complete set of toppling invariants have been explicitly constructed using the Smith normal 
form of A [15]. 

Many works have been devoted to this problem. While some of them focus on general 
properties [16, 17], others consider specific classes of graphs. When the graph is a square 
L x L grid in Z 2 , it has been shown in [15], by using toppling invariants, that G has exactly 
L non-trivial factors in (2.10). For L = 2, 3, 4, 5 for instance, one finds the following groups 

C?2x2 — Z 2 4 X Zg, G,3 X 3 = Z 2 24 x Z 112 X Z 4 , (2-H) 

G4 X 4 = Z 660 o x Z 132 o x Z 8 x Z 8 , (2-12) 
^5x5 = Z 10 2960 x Z 10 2960 x Z48 X Z 16 X Z4. (2-13) 

The structure of G has been completely determined for wheel graphs W n [18], complete 
graphs K n [19], graphs of dihedral groups D n [20], square cyclic graphs [21], and others, 
mainly Cartesian products of simple graphs. 

The identity configuration is easy to construct algorithmically. It is the only recurrent 
configuration in the equivalence class of the empty configuration hi = 0, or of any 
configuration such that Yl iGC af' = 1. A convenient choice is to start from Co with heights 
hi = Yljec^-hj = z t ~~ Zi - ft is n °t recurrent though since the conservative sites have zero 
height. However the sequence C n = Co © stabilizes for a finite index k, at which point 
Ck = Ct-i is recurrent and provides the group identity [22]. 

The identity configuration has been studied mainly for rectangular portions of Z 2 , where 
it shows intriguing and complicated fractal patterns, first explored in [22]. Its geometric 
structure in general remains poorly understood. Partial answers have been obtained for 
certain rectangular domains [23]. For square grids, very little is known apart from the fact 
that the identity on the (2L + 1) x (2L + 1) grid is simply related to that on the 2L x 2L 
grid [15]. The identity has been fully determined on certain oriented square grids [24]. 

§ If A is a non-singular matrix with entries in Z, the Smith normal form of A is the unique diagonal matrix 
D such that A = ADB for A, B in GL(7V, Z) and Dij = diSij for positive integers di satisfying di + \\di [14]. 
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Figure 1. Temperley's correspondence. The left figure shows a typical dimer configuration 
on an odd-by-odd rectangle with the upper right corner removed. In the central panel, the 
dimers touching the odd sublattice are kept, and replaced by arrows in the right figure. The 
arrows form a rooted spanning tree, where the root is represented by a circle. In ASM terms, 
there are only two dissipative sites. 

2. 4- Related models 

Through the relation between recurrent configurations and rooted spanning trees, explained 
above, the stationary regime of the ASM is connected to two other models, the dimer model 
and the loop-erased random walk (LERW). 

A dimer configuration (also known as a perfect matching in the mathematical literature) 
on a graph with an even number of vertices is a subset of edges such that every vertex 
belongs to exactly one edge. The dimer model then investigates the statistics of the dimer 
configurations. It is quite old as it was proposed in 1937 as a model of adsorption of diatomic 
molecules on a subtrate [25]. Since then it has been the subject of a vast literature, especially 
in recent years, see [26] for a recent review. 

The correspondence with spanning trees has been first observed by Temperley [27] for 
odd- by-odd rectangles in 1? with a corner removed. It is pictured in Figure 1. The dimers 
touching the odd sublattice (both coordinates are odd) are marked in blue, the others in red. 
Noting that the red dimers are unambigously fixed by the blue ones (or vice- versa), we may 
focus on the blue dimers. Replacing them by arrows, we obtain a spanning tree rooted at 
the removed corner, and living on the odd sublattice. 

This correspondence can be generalized to any rectangle and actually holds for any 
bipartite planar graph [28, 29]. If the dimers are uniformly weighted, the correspondence 
induces a uniform distribution on trees, which therefore connects to the uniform measure 
on the recurrent configurations. When the lattice has a cylindrical or toroidal geometry, 
spanning trees are to be complemented by cycle-rooted spanning forests [30]. 

The LERW has been proposed by Lawler [31] as a simpler version of the self-avoiding 
random walks (SARW). The sample paths are defined from a symmetric random walks 
by chronologically removing the loops as they form (and simultaneously shifting the clock 
backwards so that a LERW at time t has t steps). In contrast, in the SARW, the finite paths 
are those of the symmetric random walk which do not cross themselves. The finite paths 
in the LERW and SARW are identical but their construction are endowed with different 
probability measures. 

The relation between rooted spanning trees and LERW has been established in [32, 33]. 
The measure on the LERW sample paths connecting two points % and j coincides with the 
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measure on chemical paths of uniform random spanning trees between % and j. Wilson's 
algorithm [34] to generate spanning trees with uniform measure is based on this connection. 

Until recently, the simple-looking question asking the passage probability through a given 
site or a given edge had remained open, but substantial progress has been made by Kenyon 
and Wilson [35] who actually computed such probabilities on Z 2 (also on the honeycomb 
and triangular lattices). As we will see below, the LERW passage probabilities are directly 
related to height probabilities in the ASM, and lead to similarly new results for joint height 
probabilities. 

3. On representations 

A conformal theory, logarithmic or not, is primarily specified by the type (and number) of 
representations it contains. In the logarithmic case however, the situation is notoriously and 
dramatically complicated because the sorts of reducible indecomposable representations that 
may potentially appear can be of considerable complexity. The only purpose of this section 
is to briefly review and set the notation for the simplest chiral representations at c = — 2 that 
will be considered in the sequel, with no claim of exhaustivity (by far). Much more general 
results are reviewed elsewhere in this Special Issue [36]. We should also mention that most 
of the works cited in this section generalize in various ways the results given below to other 
values of the central charge. 

Highest weight representations form the simplest class. A highest weight representation 
is generated from a single (highest weight) vector \h) with conformal dimension h, i.e. 
satisfying L n \h) = (Lq — h)\h) = for n > 0. The full representation coincides with the 
Verma module Vh, obtained by acting with the negative Virasoro modes on \h), and is 
irreducible for generic values of h. 

The module Vh becomes reducible when h belongs to the Kac table, which, for c = —2, 
is the following set of rational numbers, 

[h r , s = ~ f - 1 : r, s < l} = {-|, |, f , f , . . .} U {0, 1, 3, 6, . . .}. (3.1) 

In particular V TyS = Vh ryS contains a singular vector at level rs which generates a submodule 
K+s,s> and which allows to define the quotient representation 

V riS = V rjS /V r+SjS . (3-2) 

The representations V rjS are still highest weight and generally reducible indecomposable, 
except those for s = 1, 2 or r = 1 and s even, which are irreducible. They are however 
not closed under fusion and require to consider new representations, first studied in [37], 
and shortly after in [38] where they were noted 7Z r ,i, r > 1. The representations lZ r; i are 
not highest weight representations: they are reducible indecomposable and carry a non- 
diagonalizable action of Lq. 

TZ r> i is made up of two representations Vi^r-i and Vi^r+i, whose highest weights differ 
by r — 1, tied together by rank 2 Jordan cells. They are pictorially represented in Figure 2. 
The states ip r and ip' r are the two lowest singular states in Vi^r-i, and only ip' r is set to zero. 
The situation is similar in the right representation Vi,2r+i with the two states p r and p' r . 
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Figure 2. Graphical presentation of the representations !Z r ,i- The two cone-like figures 
picture the representations Vi,2r- 1 (left) and Vi^r+i (right), where the single-colored shaded 
regions represent the irreducible sub-factors. The states ip' r and p' r are singular null states 
and generate the null Verma modules (in white). The numbers on the right show the levels 
at which the various states occur. For r = 1, the proper figure is obtained by removing the 
green region (and the state £ r ). 



The left representation Vi^r-i is a subrepresentation while the right representation 
Vi,2r+i is n °t: the Virasoro modes have an off-diagonal action from Vi^r+i to Vi^r-i- The 
main defining relations read 

L l/j r = hi^r+li'r + fr, L\ VV = 6r, L n lf) r = for Tl > 2, (3.3) 
(I^ 1 + ...)&• = Pr<Pr- (3.4) 

The pair (ip r ,ip r ) is the lowest-lying Jordan cell in TZ r ^, and gives rise to an infinite number 
of descendant rank 2 Jordan blocks lying at higher levels. 

An important feature of these representations is that the presence of the parameter (3 r . 
Once ip r is given, the fields (or the corresponding states) y? r ,£ r with their normalization are 
fixed. However ip r , as a singular descendant of £ r at level r — 1, must be equal to some 
multiple, say of a canonical choice like (Z7^ + . . .) £ r . The parameter (3 r cannot be 
absorbed in the normalizations of the various fields, and so its value is intrinsically related 
to the conformal transformation of ip r . It is therefore a number that labels the equivalence 
classes: two representations of type lZ r ,i are inequivalent if and only if they have different 
values of (3 r [39] (though all representations for fixed r have identical character). 

The parameter j3 r is known in the literature as a logarithmic coupling, an 
indecomposability parameter or a beta-invariant. A general mathematical study of chiral 
rank 2 staggered modules, of which the TZ r ,i are specific examples, has been carried out in 
[39]. The problem of non-chiral logarithmic representation is even more complicated and 
remains widely open, see [40, 41, 42] for first attempts. 

The previous scheme becomes degenerate for r = 1, because there is no field £i below 
(pi and The defining relations become (hig = 0) 

L Q ipi = ipi, L n ip r = for n > 1. (3.5) 

The graphical representation of TZ\ y i is still given in Figure 2 where the green region is to be 
removed. As a consequence, the parameter Pi no longer makes sense. 
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An important issue concerns the fusions of these two sets of representations, namely the 
quotient V TtS and the 1Z r ,i- It has been argued in [38] that the set of irreducibles V rjS (s < 2) 
and the lZ r ^ for specific values of f3 r is closed under fusion, with conjectured fusion rules 
given by 

V ril i *V ra ,i = V m ,i, V ri) i *V„ )2 = V m , 2 , (3.6) 

m m 

V ri) 2*V ra)2 = 7^,1, V ri) i* 7^,1 = 7^,1, (3.7) 

m m 

V ri ,2*K r2A = [V m -i, 2 ©2V m , 2 ©V m+ i, 2 ], (3.8) 

m 

^ri,l*^ra,l = t^m-1,1 © 2 ft m;1 © T^+i.i] , (3.9) 

m 

where the summations are \r\ — r 2 | + 1 < m < ri + r 2 — 1 by steps of 2, and Vo,2 = 7£o,i = 0. 
One sees that all T^i appear in the fusion of irreducible representations and must therefore 
carry fixed values of j3 r . The sequence of actual values starts with /3 r = —1, —18, —2700, . . . 
for r > 2, where the general term is conjecturally given by f3 r = — (r — l)[(2r — 3)!] 2 /4 r ~ 2 [43]. 
These fusions have since been confirmed on the lattice in the dense polymer model, either 
in its loop formulation [44] or in its (inequivalent) spin chain formulation [45]. The fusions 
with the other, non-irreducible representations V r , s , apart from the Vi j2r+ i, see below, are 
not known. 

Let us mention that the loop formulation of the dense polymer model gives rise to closely 
related representations, so-called Kac representations and noted (r, s), for Kac labels r, s > 1 
[46]. For each pair r, s, the characters of V r , s and (r, s) are equal, but the two are nonetheless 
distinct, except in the following cases, 

(r, 1) = V r ,i, (1, 2r) = (r, 2) = V lj2r = V r , 2 , (h.w. irreducible) (3.10) 
(l,2r + 1) = Vi j2r+ i. (h.w. reducible indecomposable) (3-11) 

For s even, (r, s) is fully reducible and reduces to a direct sum of irreducibles, but in all 
other cases, (r, s) is non highest weight, reducible and indecomposable of rank 1 (L is 
diagonalizable) . For the special value s = 3, the representation (r, 3) for r > 2 is in fact 
isomorphic to the quotient TZ rj i/ip r by the highest weight representation generated by tp r (in 
the quotient, the diagonalizability of L Q is recovered since ip r and its descendants lose their 
partner, so that the Jordan cells all become rank 1). 

All mixed fusions between Kac representations and the TZ r ,i for the special values of f3 r 
given above have been determined (conjectured) in [46]. They allow to complete the fusions 
(3.6)-(3.9) with 

ri+r-2 — 1 

n m>1 , (3.12) 

m=\ri— r2|+l,stcp 2 

V 1)2n+ i*V r2! i = (r 2 ,2n + l), (3.13) 

v — *{£:}- ,®, (3i4) 

m= \r 2 - ri --\ + - 
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Note the second line which closes on a single Kac representation. The set of representations 
V r ,2, Vi^r+i, 7tr,i) is therefore not closed under fusion. 

4. Lattice calculations 

From the general properties of the ASM reviewed in Section 2, and in view of its relations 
to spanning trees and the loop-erased random walk, one may well be tempted to suspect 
that the sandpile model is critical, in two and higher dimensions, and therefore presumably 
conformally invariant. 

The criticality of the ASM appears to related to two essential ingredients, the nature 
of the recurrent configurations and the local conservation of sand in the model (except at 
a few dissipative sites). The recurrence property of the configurations is shaped by the 
many topplings occurring in the system while approaching the stationary regime. Large 
scale topplings tend to organize far away sites in a specific way and therefore introduce 
large distance correlations. This mechanism however requires that large avalanches do occur, 
and this is only possible if most sites are conservative, enabling the excess of sand to be 
transported over large distances rather than being quickly evacuated to the sink through 
dissipative sites. It has indeed been shown that when all sites are dissipative (a uniform 
non-zero density is enough), the average avalanche size per added particle is finite and the 
correlations are exponential [47, 48, 49, 50]. On the contrary, when the density of dissipative 
sites is zero, the average avalanche size diverges with the system size like N 2 / d in d dimensions 
[7]- 

In this section, we will review the calculation of lattice correlations supporting the view 
that the ASM in two dimensions is critical and conformally invariant. Moreover we will 
collect enough evidence to argue that the underlying conformal description in the scaling 
limit is actually logarithmic. 

Before the correlations themselves, we will start by computing the marginal distributions 
of a few variables, mainly related to height variables. Although they are generally irrelevant 
for a field theoretic interpretation, it is instructive to see how they can be computed and in 
any case, they are useful to understand the calculation of the correlations themselves. 

Unless explicitly stated, all the calculations are made on Z 2 , mostly on rectangles and 
cylinders. We assume that the bulk sites are always conservative (except for a finite number of 
them when we consider the insertion of extra dissipation), whereas the boundary sites can be 
either conservative (z* — z i — 3) or dissipative (z* = ^ + 1 = 4). In the former case, the site is 
called closed, and open in the latter; we will accordingly talk about closed and open boundary 
conditions. Then the toppling matrix is given by the Laplacian, plus appropriate boundary 
conditions, Neumann or Dirichlet for closed and open boundary conditions respectively. As 
we are primarily interested in the thermodynamic limit (followed by the scaling limit), all 
results will be given in the infinite volume limit. 
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4-1. Free energies 

Perhaps the most elementary result concerns the number of recurrent configurations, which 
provides a natural definition of the partition function, 

Z = detA. (4.1) 

At finite volume (i.e. on a finite graph), Z is an integer. 

On alxM rectangle, with the two opposite edges of length L open and the other two 
closed, the partition function is given by 

L-l M 

Z i^Z} = II II t 4 - 2 cos T - 2 cos (4.2) 

m=0 n=l 

In the limit L, M — >■ 00 with r = 2^ fixed, the asymptotic behaviour of the free energy reads 
[51] 

logZr 2o Pcn , = L(M + 1) (L + M + l) log(l + V2) + \ log M + log 2f rj(r) + . . .(4.3) 

2 closed J 7f 

up to terms which vanish in the limit L, M — > oo (77 is the Dedekind eta- function). 

It follows from (4.3) that the free energy per bulk site is equal to — , and that the free 
energy per open resp. closed boundary site is ^ — \ log (1 + \/2) and ^ — \ log (1 + y/2). 
Exponentiating these numbers, we find, within the set of recurrent configurations, an effective 
number of degrees of freedom equal to 3.21, 3.70 and 2.07 for a bulk site, an open boundary 
site and a closed boundary site, to be compared with the values 4, 4 and 3 in the set of stable 
configurations. 

The result (4.3) also implies 

hm i-logZ = —L- log(l + V2) --£- + .. . (4.4) 

The third term can be compared with the general result from conformal theory. It yields 
the first and simplest indication that the conformal description has c = — 2 [9]. 



4-2. Local probabilities 

Beyond the free energy, a basic quantitative information is the distribution of the local 
values of the height variables, that is, the probability that a certain small cluster of fixed 
heights appears. On a finite lattice, these probabilities depend on the exact location of the 
cluster but become translation invariant in the thermodynamic limit. They retain however 
a dependence on whether the cluster is far or close to a boundary, and in the latter case, 
whether the boundary is open or closed. 

The simplest distribution is the probability P(/ij = a) that a given site % has height a. 
Despite the apparent simplicity, these numbers are not so easily computed. To understand it, 
let us focus on the case when the reference site % is in the bulk, far from the boundaries. The 
height can then take four values, from 1 to 4, so there are only four numbers to be computed 
(three really). The first question is: how can we characterize the recurrent configurations 
that have a given height at a given site ? 
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The answer is rather impractical in terms of forbidden subconfigurations, and becomes 
much clearer in terms of spanning trees. To see this, we use the burning algorithm in a slightly 
different (but equivalent) way, by proceeding in two steps. We first start the algorithm as 
explained in Section 2, iteratively burning all burnable sites except the site % at which we want 
to compute the height probability. When this is done, the recurrent configuration is only 
partially burnt, with a subset Cb of burnt sites and a complementary subset C u of unburnt 
sites, which contains the site % and possibly other sites. At this stage, % must be burnable -it 
is even the only one in C u to be so- and will eventually propagate the fire to the whole of C u 
(this is the second step). 

According to the burning algorithm, i being burnable means that its height is strictly 
larger than the number of its unburnt nearest neighbours (those in C u ). If this number is 
equal to k, between and 3, we have hi > k and we say that i has k predecessors || among 
its nearest neighbours. Let Xk(i) be the fraction of spanning trees such that the site i has k 
predecessors among its nearest neighbours. 

The fraction Xk(i) contributes equally to the values of F(hi = a) for k < a < 4 because 
in a recurrent configuration with hi = 4 whose spanning tree is in X k (i), one can decrease hi 
from 4 down to k + 1 and still keep it recurrent. Therefore one finds [52] 

p <". = «> = Efrp <«) 

k=0 

The number Xq(i) is the easiest and can be computed in the following way [53]. The 
spanning trees counted by X (i) have a single connection between % and one of its four 
neighbours, which can be the N, W, S or E neighbour {% is a leaf). By rotation invariance, 
which one it is does not matter in the thermodynamic limit, so we may focus on those trees 
in which i is connected to its left neighbour say, and include a factor 4. The corresponding 
spanning trees can be viewed as all spanning trees on a modified lattice, obtained from the 
original one by removing the edges from i to its N, W and S neighbours. In the modified 
lattice, the site % has degree 1, while its three N, W and S neighbours have degree 3, so that 
the toppling matrix is accordingly modified by a rank-4 matrix B as 

A = A + S[i], B[i\=(\ - 1 \ J ] on{i,N,W,S}, (4.6) 

\ 1 o o -1 / 

where B is identically zero elsewhere. By Kirchhoff 's theorem, det A gives the number of 
spanning trees in which i is a leaf and is connected to its left neighbour, while the ratio 
(det A)/(det A) is their fraction in the full set of spanning trees. It follows from (4.5) that 

1 Het A 

m = 1) = Z X (i) = ^ = det(I + A-^W). (4.7) 

Because B has rank 4, the previous determinant reduces to its four-dimensional 
restriction to {i,N,W,S}. In the infinite volume limit, A becomes the discrete Laplacian 

|| A site j is a predecessor of i if the path going from j to the root passes through i, so that j is traversed 
first. By definition, the sites of C u are the only predecessors of i. 
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on Z 2 while P(/ij = a) — >■ P a converges to a translation invariant quantity. The relevant 
entries of the Green matrix A" 1 are known exactly, from which one obtains [53] 
2(ir - 2) 

Pi = v „ ; ~ 0.07363. (4.8) 

The other three numbers Xk(i) for k > 1 are not so easy. For spanning trees in X\(i) for 
example, only one neighbour of % is a predecessor of %. So the flreline could come from the E 
neighbour of i, pass through % and go down to its S neighbour, which would therefore be the 
predecessor of %. The two bonds from % to its N and W neighbours can then not be part of 
the spanning trees in this class, and could be removed, like what we did for X (i). However 
this is no longer sufficient because the flreline coming out from the S neighbour could wander 
in the lattice and eventually touch the W neighbour of i, in which case % would have two 
predecessors at least. 

We see that there is a clear, non-local distinction for a site to have no predecessor among 
its neighbours or to have h > 1 predecessors, because the flreline can go from % to a neighbour 
directly, by using the edge connecting them, or along a path making a possibly long tour in 
the lattice. The spanning trees in X (i) are constrained by local conditions, those in Xk>i(i) 
are subjected to non-local constraints. This simple observation makes the latter much more 
difficult to compute. In fact it took over twenty years before they could be given a simple 
closed form. 

The first computation of X k (i) in the infinite volume limit was carried out by Priezzhev 
[52] by purely graph-theoretic arguments. The result however was given in terms of two 
multiple integrals which could not be evaluated analytically. The computation was later 
reconsidered in [54], which established an exact linear relation between the two integrals. 
Based on a high-precision numerical evaluation of the remaining integral, the following 
expressions for 1-site probabilities were conjectured in [54], 

^2 = --^-- 2 +- z ^ 0.1739, (4.9) 

4 Z7T 7T 2 7T d 

P 3 = - + - - ^ ~ 0.30629, (4.10) 

8 7T 7T J 

P 4 = |- ^ + A + A- 0.44617. (4.11) 
8 2n it z it A 

Remarkably, the mean value of the height distribution at a site turns out to be 
4 2^ 

(/^> = $>P a = -, (4.12) 

a=l 

a deceptively simple number in regard of the heavy technology used to obtain it, but definitely 
a hint that a simpler solution is lurking away. 

Three independent proofs for the above three probabilities were eventually given within 
a year. The first one [55] is based on a specific relation, noticed in [56], between the mean 
height (hi) and the probability that a LERW eventually visits a fixed nearest neighbour 
of its starting point, the latter being then computed in terms of dimer arrangements and 
found to be 5/16; the second one [35] also uses, in more general terms, the relation of height 
probabilities to LERW passage probabilities (see also [57] for a discussion of closely related 
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quantities in spanning trees and spanning unicycles); finally the third one [58] relies on a 
clever and direct evaluation of the multiple integrals alluded to above. The reference [35] 
especially presents a much simpler scheme to compute the probabilities P a (and more, see 
below), which provides an explanation for their simple form and the fact that they belong to 

Q[ * ] - 

Beyond the height distribution at a single site, one can consider somewhat larger clusters, 
for instance the probability P({§X2}) to have two adjacent heights 2. Generalizing the situation 
with the heights at a single site (1-site clusters), the larger clusters fall in two categories, 
those which can be computed in an elementary way, those which cannot. 

The clusters which are minimal subconfigurations (aka weakly allowed subconfigurations) 
form the easy class: a given cluster of heights is called minimal if it becomes a forbidden 
subconfiguration (FSC) when one decreases the height of any site in the cluster [53, 49]. The 
smallest minimal subconfiguration is a single height 1, whose probability has been computed 
above, while larger examples include 

<2Hi> <2H2Xi> <2X^ <|h|xi> <2hSh^ > (4.13) 

The probability of occurrence of minimal subconfigurations can be computed in exactly 

the same way as the probability to find a height 1, by cutting off appropriate bonds and 

adjusting the degrees of the sites surrounding the cluster. They are thus all given by a 

determinant of the type (4.7) where the form and the size of the matrix B depends on the 

cluster considered. In the infinite volume limit, they are given by the determinant of a finite 

matrix involving entries of the Green matrix on Z 2 (as such, they are in Q[^]). For instance, 

the probability to find the fourth cluster in (4.13) is given by [49] 

1025 1291 10549 151328 353192 429056 634880 n nnMnn , , , lS 

+ — — -= — ~ 0.000173. (4.14) 



128 8tt 8tt 2 27tt 3 27tt 4 27tt 5 81tt 6 
The height clusters which are not minimal are in general much more difficult. Presently, 

the only practical method to compute the occurrence probabilities of these clusters -and 
their correlations on short distances- is that of Kenyon and Wilson [35] , already mentioned 
above and based on the calculation of LERW passage probabilities (in addition the method 
works equally well on other lattices). Even so, the combinatorial complexity of the method 
increases very rapidly with the number of sites in the cluster, or with the distance separating 
the clusters in case of correlations. Using it, Wilson was able to compute the probability 
P(<5jK^K|H$>-) to find the heights a, b, c, d on four aligned adjacent sites and for any a, b, c, d, 
for instance [59] 

Pf/9V/qvY5\/9^ — 7042901 _ 322336899 , 28961498069 _ 30786770503 , 1294343142253 
IT^H^M^H^ - 256 512w -+- 460g7r 2 8647r 3 I0368vr4 

_ 711072840871 , 2817944526383 _ 2172951289219 , 3564889622 , 52154353024 (a i r\ 

2592tt 5 ^ 7776tt 6 8748tt 7 1 817T 5 1 2187vr 9 ' \^- lo J 

The variables we have discussed so far are related to heights. The description of the 
recurrent configurations in terms of rooted spanning trees offers other types of random 
variables. Let us recall that the oriented spanning trees can also be viewed as a loopless 
configurations of arrows, in which an outgoing arrow is attached to every site. With 
respect to the uniform distribution on all spanning trees, the probability that a certain 
local subconfiguration of arrows occurs can be computed. 
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To count the spanning trees with an arrow from i to its neighbour j, we assign the bond 
a weight 5 and take 5 to infinity to give the other bonds incident to i a relative weight 
equal to zeroP. This can be implemented by defining a new (asymmetric) toppling matrix 
A, which coincides with A everywhere except for two entries: A$j = —5 (instead of —1) and 
A^j = 3 + 6 (instead of 4). By Kirchhoff's theorem, lim^oo |det A is equal to the number 
of oriented spanning trees with an arrow going from % to j. The new toppling matrix can be 
written as a finite rank perturbation of A, as A = A + B[i,j] where the defect matrix is zero 
everywhere except on the sites where it reduces to the 2-block ' 



o o 

Likewise the determinant of the perturbed matrix A + B[ii,ji] + B[i 2 ,j 2 ] + • • • computes 
the number of spanning trees containing arrows along the edges (ik,jk)- The probability to 
have n arrows is then 

1 n 

P[arrows on (i k ,j k )\ = lim — det (l + A-'Vs^jA (4.16) 

fc=l 

It reduces to the same sort of computations as for the minimal height subconfigurations, 
although fixing arrows between neighbouring sites does not univoquely fix the heights at 
those sites. For instance the probability, in the infinite volume limit, to have two right 
arrows on the edges (i, i + ei) and (i + ei, i + 2ei) is 

P(o^o^o) = - - — ~ 0.0908. (4.17) 

It shows that the configuration —>■—>■ is slightly more likely than — >~[ or — each of which 
has probability ~ 0.0796. 

Finally let us mention that the probabilities discussed above are generally easier to 
compute on a boundary. In particular, the height distribution at a site on an open and on a 
closed boundary can be found in [60, 61]. 



4-3. Dissipation 

Closely related to the observables discussed in the previous section is the local insertion of 
dissipation. Because of the crucial role it plays in the understanding of the lattice correlations 
within a conformal picture, we give it a separate treatment. As we will see in a moment, it 
also yields the first explicit occurrence of logarithms. 

Let us consider a finite grid in Z 2 , with conservative bulk sites, the necessary dissipation 
being located on some part of the boundaries. We know that det A, with A the discrete 
Laplacian, counts the number of recurrent configuration and manifestly goes to infinity in 
the thermodynamic limit. We now ask the following question [62]: what happens if we make 
a finite number of bulk sites dissipative ? 

The dissipation rate at a site i is controlled by the diagonal entry A i;i = z* of the toppling 
matrix. A value z* — z i (— 4 for a bulk site) means i is conservative, a larger value z* > z i 
makes it dissipative. Thus a bulk site i 1 which is conservative in the original system becomes 
dissipative upon changing the diagonal entry by A iltil — > A iltil + 1, or A — > A + with 
(-Djjij = 5^5^ (in fact we could shift A iltil by any larger integer without any noticeable 

P We could assign a zero weight to the edges incident to i other than but this is less economical. 
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difference). Likewise the new toppling matrix A = A + D ix + D i2 + . . . modifies the original 
model by making the sites dissipative. 

We expect that the insertion of dissipation at isolated sites, even in finite number, 
induces non-trivial changes over macroscopic scales. When the bulk is conservative, the 
spanning trees have an arrow flow oriented towards the boundary since the boundary sites 
are the only ones to be connected to the root. That some of the bulk sites become dissipative 
-connected to the root- implies that the arrows can also flow to the root through them, 
which appear as sinks placed in the middle of the grid. The global pattern of arrows, or the 
structure of the spanning trees, is then dramatically different. 

A simple way to figure out the effect of introducing local dissipation is to see how the 
number of recurrent configurations changes, by computing the ratio det A/ det A, first at 
finite volume then in the thermodynamic limit. 

Let us begin by introducing dissipation at a single site i, far from the boundaries. The 
perturbing matrix D being of rank 1, the ratio of determinants reduces to 

^ = det(I + A- 1 A) = l + Ar/. (4.18) 

The diagonal entry A~* depends on the amount and location of dissipation in the unmodified 
model. It remains finite on a finite grid, but diverges in the infinite volume limit no matter 
what and where the dissipation is. This is to be expected because in the modified model, the 
height at site i can take the extra value hi = 5, and this allows a large number of additional 
recurrent configurations on the grid. 

The same problem persists when we introduce dissipation at several sites. One possible 
solution is to compare the number of recurrent configurations with n dissipative sites with 
the situation when there is one (rather than no) dissipative site. It is also natural to do so. 
The dissipation, if located only on the boundaries, becomes ineffective in the infinite volume 
limit because the boundaries are sent off to infinity. Thus one has to introduce dissipation 
by hand to make the model well-defined. Let us note that on the upper-half plane, the same 
remarks would apply if the boundary is closed, but not if the boundary is open. 

Therefore we compare the partition functions with n and with one dissipative site, and 
we compute det A n / det Ai. These ratios make sense in the limit because the divergences of 
the numerator and the denominator are equal and proportional to A^q. Inserting dissipation 
at n sites i±, i 2 , ■ ■ ■ leads to the ratios, 

detA ra = det ( I + A-1 ELi A fc ) = det(I + A-% g{n ,..., M 
det A x " det(I + A-iAj (1 + ArJ. J 1 ' } 

In the limit of infinite volume, this expression is actually symmetrical in the i^s since the 
denominator is translation invariant and does not depend on the position of i\. 

For n = 2, the ratio reduces to an easy 2-by-2 determinant, and yields for a large distance 
r = \ii — i 2 \ between the two insertions, 

= I logr + 2 7o + 0(r~ 2 ), (4.20) 
det Ai 7T 

where 7o = 7^(7 + § log 2) + 1, and 7 = 0.577216 is the Euler constant. 
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These first results suggest a very neat conformal interpretation. Let us assume that in 
the scaling limit, the introduction of local dissipation at a bulk site z is implemented by the 
insertion of a field ui{z,z). The ratio (4.19) would then converge to the n-point correlator 
(u)(ii) . . . u{i n )). Our prescription to normalize by det Ai and the result for n = 2 imply 

{u(z,z)) = l, (4.21) 

(u)(z, z) u)(w, w)) = — log \z — w\ + 27 , (4.22) 

7T 

while the value for n — 0, that is, no insertion at all, vanishes since it is the inverse of the 
diverging expression in (4.18), and leads to (I) = 0. 

The natural conclusion is that u is a zero weight field, logarithmic partner of the identity 
I. The 2-point function fixes the conformal transformations of uj(z, z) to be 



1 



47T 



dw 



dz 



2 



1 



u(z,z) — > u(w, w) — — log — — , or L uj — L u — — — I. (4.23) 



4tt 



Many more checks confirm this assignment [62]. Similar calculations have been carried 
out on the upper-half plane for an open and a closed boundary, and in the latter case, 
dissipation on the closed boundary + has been considered, leading to the same conclusion for 
a chiral boundary field u (x). In addition various fusions have been considered: the bulk 
fusion of oo(z, z) with itself, the corresponding boundary fusion of U0b(x) and the way u(z, z) 
expands on U0b(x) when it comes close to the closed boundary. Finally higher correlations 
with change of boundary conditions have also been computed. All of them are consistent 
with conformal correlators. 



4-4- Correlations in the bulk 

Long distance correlations are important to probe the critical properties of the model and to 
adjust its conformal description, but also raises new computational challenges. In this section, 
we review the known correlations of lattice variables located far from the boundaries, on the 
square lattice 1? . We mainly focus on the correlations of height variables, which are the most 
interesting, and briefly discuss the correlations of minimal subconfigurations and of arrow 
variables. We denote the general joint probabilities by P &...(ii, ii, • • •) = P[M«i) = °> M^) = 
b, . . .], but for the purpose of comparing with field-theoretic expressions, it is more convenient 
to use the subtracted random variables 

h a (i) = S h{i) , a - F^), a = 1,2, 3, 4, (4.24) 

and the corresponding subtracted correlations, 

P^...(M2,...) = K[h a (i 1 )h b (i 2 )...]. (4.25) 

We know that the height 1 is much easier to handle than the other three heights, so it 
is natural to start with the joint probabilities to have a height 1 at different sites. We have 
seen in a previous section that one can count the configurations with a height 1 at site % by 
perturbing the toppling matrix by a 4-dimensional defect matrix B[i], A — >■ A = A + B[i], 

+ An open boundary is dissipative in itself, so that adding extra dissipation does not change anything. Thus 
on an open boundary, the insertion of dissipation corresponds to the identity. 
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and computing its determinant. If we want to have a height 1 at different sites, we simply 
add as many blocks B[i k ] as required. Therefore, the general height 1 joint probability is 
equal to the following finite determinant, 

Pn...(ii, i2, ■■■)= det (l + A" 1 B [**]) ■ ( 4 - 26 ) 

k 

The determinant is 4n-dimensional for an ra-site joint probability, but can be reduced to 3n 
by addition and subtraction of rows and columns. When the separation distances fhi are 
large, one may expand the determinant in inverse powers of r k i = \fki\ by using the large 
distance expansion of the Green function on Z 2 [63] . 

One finds that the 2-site correlation of two heights 1 decays algebraically 

rr?n _ n 4( 7 r-2)[l + ( 7 r-2)cos4<d ( 

where r 12 = re ¥ . Only the first dominant term is rotationally invariant which suggests that 
the subtracted height 1 random variable h\(i) converges in the scaling limit to a scalar field 
with scaling dimension 2. In d dimensions, the same 2-point correlation decays like r~ 2d [53]. 
In complex notation z^i = rn, the 3-site correlation of three heights 1 is 

^111^1,^2, 13) - 5 S 7777 ~U + -.3 .3 ^2 + 



7T 9 |^2 |^i2 ^13 1 4 ^12^13^23 ^ Z 12 Z \i Z 22, 



+ 9 3( 2 7r - 2 2 \ + perm. + c.c. } + ..., (4.28) 

Z Z 12 Z \Z Z 23 J 

where the permutations to be added are the transpositions (1 -H- 2) and (1 -H- 3). The 4-site 
correlation is found to be given by, 

™ / X F t f 1 1 1 1 

Pllll(«l,«2,«3,«4) = -3- 



8 lkl2^34| 4 1^13 ^24 1 4 \ Z U Z 23\ 4 {. z \2 z Zk Z \Z z 2k) 2 
1 1 1 

+ c.c. } + ... (4.29) 



(^13 ^24 z 14 Z 2Z) 2 (^14 ^23 ^12 ^34) 2 

As we will see later, there is no need to compute higher correlations to get a convincing 
identification of the height 1 variable in the scaling limit. 

Correlations involving heights not equal to 1 are much more complicated. To date, the 
only correlations that have been computed are Pi («i,i2) for a = 2,3,4, and require heavy 
graph-theoretical and analytical calculations*. Two-site correlations for two heights larger 
than 1 still remain well out of reach. We only quote here the final results for ¥ la (i 1 ,i 2 ); the 
interested reader is referred to [64, 65] for the details. 

The large distance behaviour of the subtracted lattice 2-site correlation is found to have 
the general form 

p , iJr) = ^B+M + ^ il+ 1£a£4 _ (430) 

This situation is very likely to change. Thanks to the technique developed in [35], it is now possible to 
compute Pi a in a much more efficient way, on Z 2 and on other lattices, and could even allow to compute 
higher correlations of the type P a i...i for a > 1. 
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Figure 3. Correlations P| a (r) as a function of r, for a = 1 (blue), a — 2 (purple), a = 3 
(gold) and a = 4 (green). 



up to lower order terms 0(r~ 5 \og k r). The constants A, B,ot a , f3 a , which depend on the 
normalization of the lattice height variables and on the correlations themselves, can be 
computed analytically and take the following values 



.112 



A = -S = -2^, B = -% ( 7 + | log 2) - ^ 2 ^ , (4.31) 
ot x = 0, a 2 = 1, a 3 = 4f=2)' " 4 = _ 2§=l)' ( 4 - 32 ) 

fl _ 1 fl _ n fl — 7r 3 -57r 2 +127r-48 o _ 32+47r+7r 2 -7r 3 (a qq\ 

Pi — 1, P2 — U, P 3 — 4(7r _ 2)2 , P4 — 4(7r _ 2 )2 • (4.ddJ 

Plugging these values into the general form (4.30) yields the plots shown in Figure 3. One 
sees in particular that the heights a = 1,2,3 are anticorrelated with the height 1, implying 
that the heights 4 and 1 are positively correlated (since 2~^a^ia( r ) = 0)- 

The previous results strongly suggest that the height 1 variable converges in the scaling 
limit to a primary field of dimension 2, while the other three height variables converge to a 
logarithmic partner of the same dimension. Additional calculations will confirm this picture. 

Let us now briefly discuss the correlations of minimal subconfigurations, like those 
pictured in (4.13). Their calculation follows closely that for heights 1. In particular the 
formula (4.26) can again be used provided the various defect matrices are chosen according 
to the cluster variables one considers. In all cases, the correlation is expressed as a finite 
determinant, whose size increases with the number and the size of the clusters involved. 
These correlations have been extensively studied in [49] for 14 different clusters and their 
rotated versions. 

The correlation for any pair of minimal subconfigurations S, S', for large separation, 
takes the general form 

aa' 616; - b 2 b' 2 / 1 1 
A V? + i 

where the triplets of constants (a, bi, b 2 ) and (a', b[, b' 2 ) are characteristic of the clusters S, S'. 
The dominant term given in the previous equation is not rotation invariant, which is not 
surprising in view of the fact that the clusters themselves are not all invariant. 

The correlations of minimal subconfigurations have also been computed in the dissipative 
perturbation of the ASM, with the advantage that they contain many more terms, allowing 
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for a more reliable identification with fields in the scaling limit. A brief account of this is 
given below in Section 6.1. 

Finally the correlations between arrows follow the same calculational scheme as the 
minimal subconfigurations, as explained in Section 4.2, and are given by the similar expression 
(4.16). They are even easier to compute as the defect matrices are two-dimensional. The 
results are however different. For instance the (subtracted) correlations of two arrows, 
horizontal or vertical, separated by a large distance are given by the following expressions 

F -< 2 > = I^(7 + i) 2 + c, <^ < 435 » 

^=-i^(H) 2+o(2 ~ 4) - (436) 

^M = Teb(?-?)+^- 3 >- (4.37) 

The second correlation is obtained from the first one by a rotation by | (z — > —\z). Higher 
correlations of single arrows or correlations of larger local arrangements of arrows can be 
easily handled as well. 



4-5. Boundaries 

Boundary conditions and their changes along a boundary, are important to see how a 
statistical system responds and is an efficient way to probe the content of its conformal 
description. We have so far mentioned two boundary conditions, open and closed. We 
start by examining them before introducing and discussing two more. We will restrict to 
the simplest geometric setting, namely the discrete upper-half plane Z x N; we choose the 
boundary to be the real axis. 

The difference between open and closed boundary conditions is in the toppling matrix: 
an open boundary contains open sites ((A op )i,i = Zi + 1 = 4), a closed boundary contains 
closed sites ((A^)^ = = 3). As a consequence, A op is the Laplacian with Dirichlet 
boundary condition, A c i is the Laplacian with Neumann boundary conditions. 

It follows that one may "close" a number of sites in an open boundary by simply shifting 
the diagonal entries (A op )j ; j — > (A op )i :i — 1 at those sites. The effect of closing a segment of 
an open boundary (and vice- versa) has been studied in [51]. 

Thus we consider the ASM on the discrete upper-half plane with open boundary 
condition, characterized by a toppling matrix A op . We close a segment / of the boundary 
containing n adjacent sites; the corresponding model with the new boundary condition is 
defined by the toppling matrix A op (n) given by 

A op (n) = A op - B, for B hj = 5 iJeI . (4.38) 

The defect matrix B is the identity on the sites of / and zero elsewhere. 

The two partitions functions Z op and Z op (n), equal to the determinant of A op and A op (n) 
respectively, count the recurrent configurations for the two boundary conditions and are both 
infinite. The effect of closing the n sites of / can be measured by their ratio, which can be 
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interpreted as the expectation value of the product of two operators, (p op ' cl and (p cl ' op , flipping 
the boundary condition from open to closed and from closed back to open, 

= E [^ P ' cl (0) ^° p (n)] = dG l^ = d6t(I " A - B) - (439) 

This quantity, expected to be finite for finite n, should converge in the scaling limit to the 
field-theoretic 2-point function ((^ op ' cl (0) y? cl ' op (n))cFT of two boundary condition changing 
fields. Its exact value for n large would therefore bring important information. 

The obvious difficulty in computing the determinant (4.39) is that its size grows with 
n. The specific form of the matrix B however makes the calculation possible. As usual, 
the determinant reduces to the site indices at which B is not zero, namely the sites in / in 
the present case, where B is simply the identity. Thus the determinant to be computed is 
det(I — A~p)j Because A op and its inverse are invariant under horizontal translations, the 
determinant has the Toeplitz form det(afc_^) for k, £ the horizontal coordinates of the sites in 
/. The enormous body of results on Toeplitz determinants can then be used. 

One finds that the entries the Fourier coefficients of the following periodic 

function, 

ai(k) — Vl - cos k • |V3 - cos k - Vl - cos (4.40) 

It has the form a±(k) = (1 — cos A;) a r(A;) where r(k) is function that is single-valued, 
smooth, nowhere vanishing nor divergent on the unit circle. The asymptotic value of such a 
determinant can be computed by using a generalization of the Szego theorem due to Widom 
[66]. The result reads [51] 

Z ° P ^ d) = det(a fc _*)i< M < n ~ const, n 1 / 4 e~^ n , n > 1. (4.41) 

The exponential factor is related to the difference of boundary free energy between an open 
and a closed boundary site, equal to ^ as we have seen in Section 4.1. This non- universal 
factor must be omitted when comparing with the field-theoretic 2-point function, which 
becomes a pure power n 1 / 4 . Thus the boundary condition changing field Lp op < cl = ^p d < op is 
expected to be a primary field a dimension — |, degenerate at level 2. 

One may check this identification by considering the situation in which one closes two 
(well separated) segments on the open boundary. In this case the ratio of partition functions 
is related to the 4-point function of ip op ' c \ which can be determined exactly by assuming its 
degeneracy at level 2 (it is a complete elliptic function). On the lattice, the corresponding 
determinant is no longer Toeplitz and its asymptotic form cannot be computed analytically. 
The numerical evaluation of this determinant however compares very well with the conformal 
correlator (a complete elliptic function) and fully supports the above field identification [51]. 

Let us also mention that the converse situation, namely we open segments in an otherwise 
closed boundary, can also be examined. Despite additional subtleties, it leads to identical 
results [51]. 

The open and boundary conditions are like free boundary conditions: they refer to the 
dissipative respectively conservative character of the boundary sites, but otherwise they allow 
the heights on the boundary to fluctuate freely (within the recurrent configurations). One 
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could naturally think of some sort of fixed boundary conditions. In this class the simplest 
ones use again the language of spanning trees or arrows (remember that in the arrow picture, 
each site is assigned an outgoing arrow). The boundary conditions we want to consider force 
the arrows going out from boundary sites to be uniformly oriented, either to the left or to the 
right [67]. It turns out that this type of windy boundary conditions has a number of unusual 
features (like carrying an intrinsic orientation). It would too long to repeat the full analysis 
here, so we restrict to some illustrative steps. 

To be concrete, let us start by inserting consecutive right arrows into a closed boundary. 
From Section 4.2, we know how to force a sequence of arrows. If we want a right arrow 
between two boundary neighbouring sites i and i + ei, we simply add to the toppling matrix 
the 2-block B[i, i + ei] = ^ S 1 _(5 + 1 j . If we impose n consecutive right arrows between the 

sites i, % + ei, % + 2ei, . . . , % + ne±, we form the matrix B(5) = X]fc=o B\i + kei, % + (k + l)ei] 
and compute the following ratio, again interpreted as the expectation value of two operators 
which implement the two changes of boundary conditions, 

ML±) = E [^-(0) r ' cl W] = lim det (I + A" 1 B(5)) . (4.42) 

Zj c \ o— >oo 

It is not difficult to see that the limit over S reduces to a finite nxn Toeplitz determinant, 
whose associated function, 

a 2 ( k) = e i(fe -" )/2 Vl - cos A; • j V3 - cos k - Vl-cosA; | , (4.43) 

differs from U\{k) used above by just a phase factor. 

The phase makes a real difference though; the asymptotic behaviour of the associated 
determinant has been given in [68]. Using this result, one finds that 

= const. n- l ' A e~™ n , n » 1, (4.44) 

where the exponential factor is again due to the difference of free energy between the two 
types of boundary conditions. This result suggests that the dimensions of the two fields 
V? c1 '^, y^' cl add up to —1/4. If we restrict to the dimensions in the Kac table, the only 
possibility is that the two dimensions are —1/8 and 3/8 (the two fields need not be equal as 
the arrows are outgoing from the closed boundary in one case, incoming in the other case). 
But with these two values, we seem to run into trouble since their conformal 2-point function 
vanishes identically. The apparent paradox can be solved if we remember the essential role 
of dissipation. 

When explicit calculations are carried out on the infinite planar lattice, one generally 
starts with a finite system and takes the infinite volume limit of the results obtained at finite 
volume. As we know, this is well-defined provided the finite system involves dissipation, 
usually located on the boundaries. Even if the size of the system is ever increasing, dissipation 
is maintained at every step. However the conformal formulation on the full plane describes 
the system which is rightaway in its infinite volume limit, and so by itself is not aware that 
disspiation was present. Thus in order to include explicitly the presence of dissipation on the 
far away boundaries, one has to insert the field w(oo) representing the insertion of dissipation, 
as discussed in Section 4.3. 
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In the present case, the boundary being fully closed, we must add dissipation by hand, 
so that the power law in (4.44) should be identified in the scaling limit with the 3-point 
function given by {oj(oo) (p cl, ~*(0) M ,cl (n)). Its general form is indeed given by the correct 
power ~ n~ 1//4 . 

The same analysis in an open boundary leads to two further boundary condition changing 
fields, ip op ^ and M ,op , both of dimension 0. 

To go further, one may look at 3-point functions. Namely we change the boundary 
conditions three times by inserting two stretches of different conditions, of lengths say n\ and 
n 2 . The ratios of partition functions can be computed as above, by using the appropriate 
defect matrices, and all take the form of a finite determinants. As these do not have the 
Toeplitz form, they must be evaluated numerically. 

If we assume that the boundary conditions changing fields are primary, the generic form 
of a 3-point function is 

(MO) M*) <Ps(v)) = -nst JlZTy^Z^ < 4 - 45 ) 

By varying the lengths x = n\ and y = rii+ri2, a fit on the numerical data allows to determine 
the conformal dimension of each field. For this purpose, numerical calculations have been 
carried out for n\ = 20, 30, 50 and 70, and for each of these values, n 2 was varied from 10 to 
150. 

Two instructive situations are pictured below. 

ni ri2 n\ ri2 



op cl op op cl op 

For these two cases, the numerical data show that the related 3-point functions}! are 
consistent with the forms ( ^^i ) 1 / 2 and n^ 4 respectively. Comparing with (4.45) determines 
or confirms the following weights, 

h op ' cl = h cl ^ = -i h op ^ = h op ^ = 0, h cl ^ = §. (4.46) 

The last cases to consider is when the direction of arrows is changed, from right to left 
or vice- versa. The two changes of orientation are in fact different. 

When right arrows are changed to left arrows, the two opposite arrows o— >■ o <=-o point 
to a same site i. Whether % is open (connected to the root) or closed makes a substantial 
difference. Recalling that the arrows eventually flow towards the root, an open site i allows 
the two sequences of arrows along the boundary pointing to i to go directly to the root 
through i. If i is closed, the flow has to go back in the bulk of the upper half-plane to find its 
way to the root. Thus the presence or absence of dissipation at the single site i has strong 
effects on the global pattern of the arrow flow making up the tree. Thus we have to consider 
two different boundary condition changing operators, M°^^~ an d ■ 

When the arrows are changed from left to right o <=— o o— > o, the previous duplication is 
not necessary since the sites concerned are not dissipative anyway (their arrow cannot point 
to the root by construction). Thus a single field is sufficient. 

jj Note that in the two cases, most of the boundary is open, i.e. dissipative, so there is no need to insert 
extra dissipation. Therefore in the conformal picture, the two situations are described by 3-point functions 
with no dissipation field u inserted. 
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«- 
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hi] G Vl >2 


[0] G Vi )3 


I 



Table 1. Representations containing the fields (p a 'P which implement a change of boundary 
condition from a (row label) to f3 (column). The numbers in square brackets denote the 
scaling dimensions. 



By following the same steps as above with 3-point functions, one obtains 

h -+f<- = h ^ = o, = 1, (4.47) 

completing the list of conformal weights for all possible fields switching among the four 
boundary conditions. 

Beside their conformal weight, a much more important issue concerns the type of 
conformal representation these eight fields belong to. A convenient way to answer it is 
to use the constraints coming from the composition law of all these fields, expressed by the 
boundary fusion algebra. Namely, the fusion <y3 a ' 7 *<^ 7,/3 must close on fields which interpolate 
between the boundary conditions a and /3, and should in particular contain the boundary 
condition changing field ip a ^ . 

Using these constraints, a proposal has been made in [67] that identifies the specific 
representations accomodating the eight fields, on the basis of two main assumptions: (i) the 
boundary condition changing fields are primary, degenerate at a level as low as possible, and 
(ii) they belong either to a highest weight representations V TtS or to a rank 2 logarithmic 
representation TZ r ^ of the type reviewed in Section 3, or to a quotient thereof. The results 
are summarized in Table 1. 

To take an example, consider the field (p° v ^ . From the two assumptions, (^ op ' cl and (p cl '~* 
of weight —1/8 belong to an irreducible representation V lj2 . Therefore ip op ^ should occur 
in the fusion Vi j2 * Vi j2 = IZi^ (see (3.7)). The representation TZi^ contains two fields of 
weight 0, a field <pi that behaves like the identity, and its logarithmic partner ip\ (see Figure 
2). The field ip op '^ cannot be the identity, hence must be identified with ipi. However the 
primary partner ipi must be null since the identity does not interpolate between two different 
boundary conditions. It follows that (p op '^ should belong to the quotient lZi,i/(pi = Vi )3 , and 
is degenerate at level 3 (it can be shown that its descendant at level 1 is not null [67]). 

Table 1 can be completed using similar arguments [67]. The representations mentioned 
are those obtained by using the fusion rules as reviewed in Section 3. Additional mixed higher 
correlators have been computed to cross-check the identifications of Table 1, namely 4-point 
correlators involving four changes of boundary conditions, or three changes of boundary 
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conditions and the insertion of dissipation on a closed boundary, or of a boundary height 1 
field. All of them have been found to be consistent with the proposal in Table 1. 

However we cannot exclude the possibility that the fields in fact belong to quotients of 
these representations. Indeed we have observed in [67] that the primary fields tp% and tp^ 
belonging respectively to 7?-2,i and decouple in some correlators without being able to 
prove that they are actually null. If this turns out to be the case, the corresponding boundary 
condition changing fields would belong to quotients T^i/v^ and TZs t i/ips which are no longer 
logarithmic, as explained in Section 3. 

On the other hand, if the fields ip 2 and y? 3 are not null, some of the boundary condition 
changing fields would belong to logarithmic representations, of rank larger than 1. The 
physical significance of this and the physical interpretation of the logarithmic fields remain 
open questions. 

4-6. Boundary effects in height probabilities 

The height probabilities at a single bulk site have been discussed in Section 4.2. At an infinite 
distance from all boundaries, they are given by the four numbers P a discussed in Section 4.2. 
Although they are probabilistically and combinatorially interesting, they are not so relevant 
for the conformal perspective. Much more important and relevant are the corresponding 
one-site probabilities at a finite but large distance to a boundary, especially if we impose 
different boundary conditions. 

In this section, we review the exact results for the 1-site height probabilities F a {m) on 
the upper-half plane with open and closed boundary conditions, where m is the distance to 
the boundary. Naturally the limit \im. m ^^¥ a {m) reproduces the numbers P a . 

When the boundary condition is homogeneous on the real axis, these probabilities reflect 
the precise nature of the fields associated to the height variables since they correspond to 
1-point functions of bulk fields, equivalently to 2-point functions of chiral fields. Moreover, 
the probabilities for open and closed boundary conditions are related to each other by the 
appropriate insertion of the boundary condition changing field (p op ' cl and therefore provides a 
highly non-trivial check of the consistency of the conformal picture obtained so far for these 



The graph theoretical techniques to compute the functions P (m) are basically identical 
to those used in the full plane, with two extra complications. First the full translation 
invariance is lost and reduces to the horizontal invariance only; this increases the number of 
graphs that need be handled separately. Second the functions P a (m) are obtained as multiple 
integrals depending on the distance m which require a rather long asymptotic analysis to 
compute the dominant contributions. 

The following results have been obtained in [69, 54] (and first in [70] for a = 1) when 
the boundary, here the real axis, is either open or closed, 



fields. 



P° p (m) = P a + \{c a + ^ + d a logm) + . . . , 
P c a '(m) = ^ - \{c a + d a logm) + . . . , 



(4.48) 



(4.49) 
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up to terms of order 0(m 3 log fc m). The coefficients are explicitly given by 



c i = T = rf i = 0, ( 4 - 5 °) 

^ = 1^(7 + flog 2)-^, d 2 = ^, (4.51) 

C 3 = ^(7 + |log2)+^^, d s = & (4.52) 

c 4 = -(ci + c 2 + c 3 ), d 4 — —(di + d 2 + d 3 ). (4.53) 
One can make two observations regarding these results. 



The first one concerns the form of these probabilities. The functions Pi (to), for 
both boundary conditions, are algebraic while all the others P a >i(m) have an additional 
logarithmic term. It follows that F 3 (m) and P 4 (to) can be written, to order to" 2 , as linear 
combinations of Pi (m) and F 2 (m). Setting 

F a (m) -P a = a a [F 2 {m) - P 2 ] + /3 a [Pi(m) - Pi], (4.54) 

we obtain the coefficients 

«i = 0, a 2 = 1, a 3 = 2§=5j, "4 = -2§J%> ( 4 - 55 ) 

/3i = l, & = 0, /3 3 = ^" 5 4 g+gr 48 , /3 4 = 32+ 4 4 ffgr 3 , (4.56) 
that is, exactly those given earlier in (4.32) and (4.33) ! We stress that the linear combinations 
(4.54) hold for either boundary condition with the same coefficients. 

Second, the fact that the same set of coefficients c a , d a control the probabilities for the 
two boundary conditions confirms that they are closely related. Precisely how depends on 
whether the open/closed condition changing field has been properly identified as a primary 
field of weight —1/8, but more importantly, on the exact nature of the fields which are 
supposed to describe the four height variables in the scaling limit. 

This should provide a crucial test to understand not only the above probabilities on 
the UHP, but also the bulk correlations reported in Section 4.4 and the striking equality of 
the coefficients a a , (3 a . We believe that the successful understanding of these specific feature 
constitutes one the most convincing support for the logarithmic conformal picture. The 
details for this are given in the next section. 

5. Conformal height variables 

The natural microscopic random variables of the Abelian sandpile model are the height 
variables h(i), assigned to the vertices and taking the four values 1,2,3,4. The distribution 
of h(i) at a single site and in the infinite volume limit can be computed exactly (see Section 
4.2) but is of little value from the conformal point of view. Since conformal invariance enforces 
zero expectation values for the fields with non-zero scale dimension, one is more interested 
in the subtracted height variables. The joint probabilities we have computed so far show 
that the variables h a (i) pertaining to each of the possible height values at site i make sense 
separately, and that they presumably converge in the continuum limit to conformal fields 

h a (z,z), 

K(i)=5 h(:i):a -F a — > h a (z,z), a = 1,2,3,4. (5.1) 
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The main issue is now to guess enough of the conformal nature of the fields h a (z, z) so as to 
reproduce all the lattice calculations reviewed in Section 4. 
From these, we learn a number of important facts: 

• all four fields have conformal weight (1,1), so a total scale dimension 2, 

• hi is a non-logarithmic field, whereas the other three h a> \ are logarithmic, 

• looking back at (4.50) and (4.51), we observe the equality d 2 — c±, namely the logarithmic 
term of P 2 (m) has the same coefficient as the algebraic term of Pi(m), a signal that 
h 2 (z,z) could be the logarithmic partner of hi(z,z), 

• the linear relation (4.54), confirmed by the bulk correlations (4.30), suggests that h 3 , h 4 
are linear combinations of hi, h 2 , namely h a (z, z) = a a h 2 (z, z) + (3 a hi(z, z). 

It is therefore sufficient to find the conformal nature of two fields, hi(z, z) and h 2 (z, z). 

The following conjecture has been first formulated in [69], and further analyzed in [54]. 
It does not completely specify, and by far, the full representation to which hi,h 2 belong, 
but rather sets the minimal framework needed to compute the relevant correlations. The 
conjecture is as follows: 

The two bulk height fields (hi,h 2 ) = {tp,ip) form a logarithmic pair of weight (1,1). The 
field tp is primary and degenerate at level 2, while the conformal transformations ofip read 

L ip = L ip = ip + \tp, Liip = p, Litp = p, L n> ii) = L n> iip = 0, (5.2) 

L n > p = L n > p = 0, L p = p, L p = p, L n >ip = L n >ip = kI 5 n: i, (5.3) 

where the two fields p and p have weights (0,1) and (1,0) respectively, and k is a computable 
parameter determined by the normalization ofip. 

The previous transformations are all we need to compute the seeked correlations. They 
are reminiscent of the relations defining the representations TZ rt i of Section 3, and can be 
seen as a partial characterization of a non-chiral extension of 7Z 2> i (the three fields ip, ip, p 
satisfy the defining relations of 1Z 2 ,i with respect to the left Virasoro modes, as do ip, tp, p 
with respect to the right modes). We will see that they are also consistent with the two 
additional relations 

L_ip = L_ 1 p = l5\p. (5.4) 

As in the chiral case, the parameter j3 ought to be an intrinsic parameter that labels 
inequivalent representations [42]. On the other hand, the parameter A is clearly related 
to the normalizations of tp and ip, which are fixed if we insist that tp and ip are the scaling 
limit of the height 1 and height 2 lattice variables. 

Altogether this yields seven fields {I, p, p, dp, dp, p, ip} with total scale dimension less 
than or equal to 2. We stress again that this does not imply that these seven fields are the 
only ones on the lowest levels. 

On the basis of the assumptions given above about the conformal nature of the fields 
ip, p, we can compute the correlators and compare with the results obtained on the lattice. 
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5.1. Height variables on the upper-half plane 

The first correlators we would like to compute are (fi op,cl (zi) fi cl,op (z 2 ) tp(z, z)) and 
fj cl ' op (z 2 ) ip(z, z)), where zi, z 2 e R lie on the boundary whereas z is in the bulk of 
the UHP. Since ip and ip are identified respectively with hi and h 2 , the height 1 and height 2 
variables, the two correlators, properly normalized by {fj op ' cl (zi) fi cl ' op (z 2 )), should correspond 
to the probability that the site at position z has height 1 or 2 when the boundary condition 
on the real axis is all open except on the interval [zi, z 2 ] where it is closed. In particular the 
two limits 

lim (/i op ' d (^)/i cl - P(^)^(z,^)) aM Um {^\ Zl )^ p {z 2 )^{z,z)) 

zi^z 2 (/i°P' cl (^i) fl d > op (z 2 )) -*i,«^oo (fl°P> cl (zi) fl cl > op (z 2 )) ' 

should reproduce respectively the functions F 2 p (m) and F^m) given in (4.48) and (4.49) 
for z — z* = 2im. Similarly the same limits with ip replaced by p ought to reproduce the 
functions P° p (m) and Ff(m). Then P° p (m) and F d (m) for a > 3 will follow from the linear 
combinations h a = a a ip + (3 a (p. 

Before going on, let us note that the two limits are in fact different. The first one 
leaves the boundary fully open, while it is fully closed in the second case. As follows from 
Section 4.5 where we have stressed the role of dissipation, the numerators should respectively 
correspond to (ip(z,z)) op and (u(oc)ip(z,z)) c \. This is not inconsistent with the fusion 
Vi,2 * Vi 5 2 = 7?.i,i if we interpret it in two different ways. In the first limit, the 
boundary fusion /! op ' cl */i cl ' op must close on fields living on an open boundary. We have noted 
before that the dissipation field on an open boundary corresponds to the identity, so that 
the fusion does not close on T^i but on the quotient Tl^i/ipi = Vi,i. The second limit deals 
with the fusion /i cl ' op * /i op ' cl leaving a closed boundary on which the dissipation field is not 
trivial, so that this one really closes on fields in TZi^. Therefore we have [62] 

^ p ' c \z)^ op (0) = z 1 / i C I ^I + ... (5.6) 

^ op {z) /i°P> cl (0) = z l " HO) + A I log*] + . . . (A = -i). (5.7) 

reinforcing the consistency of the conformal picture. 

If the calculation of (nn^p) is easy since both fi and tp are degenerate at level 2, that 
of (fJ.fJ.ip) is tedious because the conformal transformations of ip are not homogeneous (both 
are non-chiral 3-point functions but reduce to chiral 4-point functions). It turns out that the 
most general form of (fifiip) depends on 13 arbitrary coefficients, which can be reduced to 3 
by imposing appropriate physical requirementsff. Two of them are related to the norm of ip 
and to the multiple of <p which can be freely added to it; when ip is identified as the scaling 
limit of the height 2 variable, these two coefficients are related to the numbers c 2 , d 2 obtained 
from the lattice. Interestingly the third coefficient is related to (3 through the relations (5.4). 

For (3 7^ 0, one finds the following two correlators [54] 

(fi( Zl )fi(z 2 )^(z 3 , z A )) = ^f- X pr ^- hc 2 + 2d 2 log ^ + ^ - ^ X / —?- 

z 34 V 1 — x I ^ ^ 4 V 1 — x 

ff These are not related to the conformal nature of the field tp: one requires that the limit 
lim Zl ^ Z2 z 1 2^ 4 (fJ-(zi)fi(z2)'>p(z, z*)) does not contain a logarithmic term logzi2 and depends only on z — z* , 
see [54]. 
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+ d 2 



z 2 r 1 1 



^13^24 



+ 



x-2 2y/T=x 



and 



+ * — H7 7r^ 1 1 + 3 log M " 2 T3ir - ^ j • (5 ' 8) 

^ z 1 / 4 x -2 

(p(z 1 )p(z 2 )<p(z 3: ^4)) = — ^ -3 — 7^=' ( 5 - 9 ) 

A ^ 34 VI - x 

where z^- = Zi — Zj, x = fj^f^ and z± = z\. The constant k is equal to k = — d 2 . 

Let us see how the above results reproduce the probabilities computed on the lattice. 
We set z 3 = z and z 4 = z*. 

We have to look at the limits when z±, z 2 coincide, either at a finite point of the real axis, 
say the origin, or at infinity. So we set —z\ = z 2 = R and see what happens when R — > or 
+00. From 

(R-z*)(R + z)' v ; 

we see that 1 — x has norm 1, and makes a full circle around when R takes the values from 



to +00. This implies that \/l — x goes to +1 when R — >■ and to —1 when R — > +00. As 
the above correlators have this term as prefactor, one expects a change of sign in the two 
limiting cases, corresponding respectively to an open or a closed boundary. This is also the 
case in the probabilities given in (4.48) and (4.49). 

Another important consequence of the change of sign is that the very last term in (5.8) 
is singular when R — > 00, and brings a singularity log2 12 . As the limit of the correlator is 
supposed to reproduce ¥f(m), which is perfectly regular, the singularity cannot be present. 
It either forces d 2 = 0, but then all logarithmic terms log 234 ~ logm disappear, or else the 
value of the (3 parameter is 

Pasm = \- (5.11) 

With this value an easy calculation shows that the limits of z[ 2 x ^ (p(zi) p(z 2 )ip(z , z*)) in (5.5) 
exactly reproduce the two expressions (4.48) and (4.49) for a = 2. 
The same limits of the other correlator yield 

H *•)> = or (5.12) 

and again equal (4.48) and (4.49) for a = 1 provided A = — \. 

To the best of our knowledge, the ASM provides the first lattice realization of a 
representation with the value j3 — \. It is well-known [71] that the bosonic sector of 
the symplectic free fermion theory contains fields pg, pg, ipg, ipg with identical conformal 
transformations to the above p, p, tp, ip but a different value of the parameter, namely fig = — 1. 
An extensive comparison between the two cases has been made in [54]. The two values 
j3 = — 1, \ are precisely those which are discussed in the Example 7 of 
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5.2. Height variables across a strip 



The formula (5.8) and (5.9) allow to test the effect of a change of boundary condition and 
to compute finite-size corrections which can then be compared with the results of numerical 
simulations. 

A case which is particularly instructive is when the UHP is bordered by a boundary 
which is closed on the negative real axis and open on the positive part. Setting z\ = — oo 
and z 2 = in (5.8) and (5.9) and taking j3 — \ and A = — |, we obtain the conformal 
prediction for the probability that the site z in the bulk of the UHP has height a in presence 
of a half-closed and half-open boundary 

lim z[^ (/! op ' cl (^) pt cl ' op (z 2 ) K(z, z)) 



Pf*(z)-P a 



Z\ — > — OO , 22— >0 

z + z* 



|2c a + 2rf a log 



+ 



i u d a z + z 

T \z\ 



(5.13) 



\z\(z-z*) 2 

where the coefficients c a ,d a are those appearing in (4.48) and (4.49). 

The conformal map w = ^ log z transforms the UHP to an infinite strip of width L, 
under which the negative, closed, and the positive, open, parts of the real axis are mapped 
respectively onto the lines Imw = L and Imw = 0. The transformation of the (non-chiral) 
3-point function (fifih a ) then allows to transport the previous probabilities to the strip. It 
requires the finite conformal transformation law of ip, which can be computed by integrating 
the infinitesimal transformations (5.2) and (5.3). The result reads, for A = — |, 

ipuh P (z,z) - J log \z'(w)\ 2 <f uhp (z,z) 



Atvip(w,w) = \z'{w) 



+ 



z"(w) 



z"(w) 

Puhp(^, Z) + Puhp(z, z) + K 



z"(w) 



2z' 2 (w) 



(5.14) 



2z' 2 H p 2z' 2 (w) 
Using this, we obtain the following expression for the probability on a strip with 
coordinate w = u + if , < v < L, 



K - Pa = (j) ~ 

\LJ si 



7T\ 2 COs(7Tf / L) 

sin 2 (7Tf / L) 



1 + cos(— ) 

hi 



+ d a \og -sm(-) 
In L 



+ 



dqTX 2 

4L 2 



(5.15) 



If the field identifications conjectured earlier are correct, this formula gives the (subtracted) 
probability, on an infinitely long strip of width L, that a site at position w has a height 
value equal to a. By translation invariance, it only depends on the transverse coordinate v, 
which is also the distance of the site to the open edge, L — v being the distance to the closed 
boundary. 

The previous probabilities would not be easy to compute analytically on the lattice, 
but they can be estimated from numerical simulations (for a recent reviwe on numerical 
simulations in models like the ASM, see [72]). For these, the ASM dynamics is run on a 
finite rectangle of height M and the four probability profiles are measured across the strip. 
For a ratio M/L large enough, the numerical profiles should be good approximations of 
what they are on a infinite long strip, themselves presumably given by the above conformal 
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Figure 4. These plots show the subtracted probabilities to find a height 1 (orange), 2 (in 
blue), 3 (green) or 4 (red) at a site lying along a transverse line of a rectangle M x L, 
and going from an open boundary (left end of the graph) to a closed boundary (right end). 
The dots correspond to data obtained from numerical simulations, while the solid curves 
represent the conformal predictions (5.15). One notes the characteristic change of sign when 
going from an open to a closed boundary, discussed in Section 5.1. 



predictions. To obtain the numerical data, about 10 10 recurrent configurations have been 
sampled on a rectangle with M = 200 and L = 50 [69]. 

The results are given in the plots presented in Figure 4 and show a remarkable 
agreement. They provide very strong support to the conjecture identifying the four height 
variables h a (z, z) as linear combinations of the two partner fields (p, ip, having the conformal 
transformations given in (5.2), (5.3). 



5.3. Height correlations on the plane 

Whether the field assignment as well as the conformal transformations conjectured for the 
height variables can reproduce the correlations computed on the plane is another consistency 
check. It is non-trivial for the following reason. From the discussion in Section 4.5 about 
the insertion of dissipation, it follows that the proper way to understand the lattice 2-height 
correlations computed on the plane is to say that they should be given, in the scaling limit, 
by a conformal 3-point function, 

P q6 (zi 2 ) - Vb = (h a {zi, zt) h b (z 2 , z 2 ) w(oo)>, \z 12 \ > 1. (5.16) 

Note the happy by-product that the correlation of two heights 1 is then not identically zero, 
as would have been the case had the field uj not been inserted. 

When a or b is equal to 1, the probabilities P a b(^i2) are known to dominant order and 
given in (4.30) so we can check whether their forms are consistent with the 3-point conformal 
correlators. For a,b > 1, the previous equation is a conformal prediction. 

The 3-point correlator (h a (zi, z\) hb(z 2 , z 2 ) u(z 3 , z 3 )) can be computed from the 
conformal Ward identities or directly by applying the finite Mobius transformation that maps 
the three points zi,z 2 ,z% onto three fixed points, say —1,-1-1,0. From the transformations 
laws (4.23) and (5.14) of u and ip, one easily obtains 

iftzi, z x ) = + log |^| p(-l) - p(-l) - p(-l) - «}, (5.17) 
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1 



2tt 



^12 



w(z 3 , z 3 ) = w(0) - — log — ^- . (5.19) 



Inserting these in the correlator, one sees that it retains a logarithmic singularity when 
z 3 — > oo unless all 2-point functions of ip and ^ vanish identically, 

(ip(z!, z 1 )ip(z 2 , z 2 )) = ((fi(z!, z 1 )i)(z 2 , z 2 )) = (<f(z!, zi)<p(z 2 , z 2 )) = 0. (5.20) 

If the last equation is familiar, the first two are more peculiar and mark another difference 
with respect to fields ipe and ipe discussed earlier. 

One then obtains that the mixed correlators are related in the following way, 

(<p{z x , z 1 )ip(z 2 , z 2 )u(oo)) = -, (5.21) 

| Zi 2 1 

((pizuzMfa, hMoo)) = y^—ttIb + Alog \zi 2 \\ (5.22) 

\Z\ 2 \ < > 

<V(*i, ^1)^(22, Z2M00)) = t^tAc + 2B log \z 12 \ +A\og 2 \z 12 \\. (5.23) 

\Zl2\ 1 J 

Taking the appropriate linear combinations yields the required correlators [54], 
(h a {z 1 , zi) h b {z 2 , z 2 ) w(oo)) = Y^-rrlaa^bC + (a a (3 b + (3 a a b )B + f3 a f3 b A 

\Zl 2 \ < 

+ [2a a a b B + (a a f3 b + (3 a a b )A] log \z 12 \ + a a a b A\og 2 |^i 2 ||. (5.24) 
For 6=1 and from the values of the coefficients cci = 0, /?i = 1, it reduces to 

(/i a (zi,^i) h x {z 2 ,z 2 ) w(oo)} = -^--{a aJ B + f3 a A + a a A\og \z 12 \\. (5.25) 

\Zl2\ L J 

This exactly matches the form of the 2-site probabilities (4.30) computed on the lattice. The 
three constants A,B,C can only be obtained from lattice calculations; the coefficient C is 
presently unknown. 

Numerical simulations have been carried out to test the 2-site probabilities F ab (r) when 
both a, b are larger than 1. A fitted value of C ~ —0.009 shows an excellent agreement, even 
on surprisingly small distances (see [54] for plots). 

Let us recall from Section 4.4 that when they involve the height 1 variable only, more 
correlators are known. In fact it has been observed [49] that the scaling limit of the 2-, 3- 
and 4-site probabilities are all correctly reproduced if the height 1 variable is identified in the 
scaling limit with the following composite field in the symplectic free fermion theory [73, 74], 

h 1 (z,z) = <p (z,z) = -F 1 :de8e + ded6:, Pi = 2 ^~ 2 \ (5.26) 

provided the zero weight field ue —: 09:, which plays the role of dissipation, is inserted in 
the field theoretic correlator. In particular the 4-point correlation {(p g (l)(p e (2)(p g (3)ug) = 
vanishes identically, in agreement with the result that the lattice 3-site joint probability 
decays with a global power 8 in the distances (rather than a power 6), implying that its 
scaling limit vanishes. 
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More generally, it has been shown [49] that the (subtracted) variables related to minimal 
sub configurations S can all be identified with the following combination of fields of weights 
(1,1), (2,0) and (0,2), for appropriate coefficients a,bi,b 2 depending on S, 

h s (z,z) = - ja -.8989 + 0686: +b x -.8989 + 8989 : + ib 2 -.8989 - 8686:} . (5.27) 

For each of about a dozen minimal sub configurations, the three coefficients have been 
determined exactly, and all mixed 2-point correlators have been computed. They all match 
the lattice result given in (4.34). This conclusion holds in much greater generality, since it 
has been proved in [75] that the variable associated to any local bond modification converges 
in the scaling limit to a field of the form (5.27). A more direct derivation of these fields from 
the defect matrix used to compute the correlations has been given in [76]. 

Thus even though the variables associated with the minimal subconfigurations, including 
the height 1 variable, can be consistently described within the free symplectic fermion theory, 
it turns out that this is no longer the case for the higher height variables which force a different 
value for the parameter (3. 

6. Other developments 

In this last Section, we would like to briefly mention a number of complementary results 
and further developments. In most cases they complete the current picture we have of the 
Abelian sandpile model. 

6.1. The massive sandpile model 

In the lattice calculations we have presented so far, the dissipation is localized on the 
boundaries or at isolated sites in the bulk. In all these cases, the density of dissipative 
sites vanishes in the infinite volume. We have argued that this makes the system correlated 
over large distances, with the consequences that the avalanche size distribution after the 
addition of a single grain has a power-lawed tail, and the correlations decay algebraically. 

As mentioned before, the situation is drastically different when all sites, or at least a non- 
zero density, are dissipative. It is known in this case that the average size of the avalanches 
is finite and that the correlations decay exponentially [47, 48, 49, 50]. 

The dissipative model is defined simply by changing the diagonal entries of the toppling 
matrix. The simplest possibility is to choose A^j = z* = z i + t for all sites, including the 
bulk sites. As compared to conservative bulk sites, for which z* = z i , the height variables 
now take t extra values, namely they range between 1 and z i + t (so between 1 and 4 + 1 on a 
square grid), and more importantly, each time a site topples, t grains of sand exit the system. 
This causes a substantial loss of sand during the relaxation process, and thus weakens the 
avalanches which then travel over typically much smaller distances. The conservative model 
is recovered in the limit t — > 0. 

The model with t > is non-critical. The correlations are controlled by the inverse of 
the toppling matrix, a discrete massive Laplacian, whose large distance behaviour is given by 
the modified Bessel function K (ry/i) ~ e _rv/ *, implying that the correlation length diverges 
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at t — like £ ~ l/\/i. The scaling regime is reached by taking simultaneously the large 
distance limit r = - — > oo and the critical limit t = a 2 M 2 — > in terms of a scale a — > 0. 

a 

The product ry^ — > Mz then defines an effective mass and a macroscopic distance. Thus 
the insertion of dissipation everywhere in the sandpile model is like a thermal perturbation 
away from the critical point, and corresponds to a massive theory. Around the critical point, 
the average avalanche size is finite and behaves like 1/t [48]. 

The correlations of minimal subconfigurations in the dissipative model have been 
computed on the square lattice [49] . For instance, the correlator of two heights 1 is given at 
dominant order by 

PnM = ~t 2 P? - \K Q K<< + + l 4^Kl) + ... (6.1) 



2 u 2 u u ' 2vr u ' 4tt 2 

up to higher order terms in t, and where the argument of all Bessel functions is ry/t. In the 
prefactor, F 1 is the height 1 probability at the critical point. 

As the height 1 has the field identification (5.26) in terms of massless symplectic fermions, 
it is tempting to extend it to the massive case. Indeed the scaling limit of the correlator (6.1) 
exactly matches the 2-point function of 



tp 9 (z,z) = -Pi 



2 



-.8989 + 8989: + ^- :99: 



(6.2) 



where the two symplectic fermions are now massive and described by the action S = 
\ J 8989 + *g-96, with 2-point function (9(z)9(w)} = K (M\z - w\). 

Similar results have been obtained for the variables related to minimal subconfigurations 
[49]; nothing is known however of the massive extension of the higher height variables h,2, h% 
and h±. 



6.2. Boundary height variables 

We have discussed at length the height variables on the plane or on the upper half-plane, to 
find that in the bulk, the height 1 has a very different scaling behaviour compared to the 
height 2, 3 and 4 variables. A natural question concerns the nature of the height variables 
on a boundary. This has examined for open and closed boundaries, in [60, 77] for the critical 
model and in [61] for the massive model. 

That fact that the height correlations can be computed much more easily on a boundary 
than in the bulk has been noticed in [60], where 2-site correlators have been first obtained 
in the critical model. However they do not provide enough information to decide whether or 
not the different height variables have different scaling properties. This question was settled 
almost simultaneously in [77] by the calculation of multisite correlators in the critical case, 
and in [61] by looking at 2-site correlators in the massive model, which allows for much finer 
field identifications. 

The results depend on the boundary condition, open or closed, since the height variables 
take a different number of values, given in the critical case by 4 and 3 respectively. The main 
conclusion is that, for both boundary conditions, none of the height variables is logarithmic. 
Again the field identifications can be more conveniently expressed within the symplectic free 
fermion theory. We simply quote the results in the critical model. 
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For an open boundary, one finds that all four height variables are proportional to the 
same field, 

hl p (x) = N° p :d6d6 : , a = 1,2, 3, 4, (6.3) 

where the normalizations N° p are known exactly, and where the fields 9, 9 satisfy Dirichlet 
boundary conditions. 

On a closed boundary, the three height variables converge to different scaling fields, 

h*(x) = : 8989 : + M a cl : 9889 : , a = 1, 2, 3, (6.4) 

for known coefficients iV^ 1 , M^ 1 , and where the fermions satisfy Neumann boundary 
conditions. In particular Mf = so that the height 1 variable behaves the same way 
on the two types of boundaries. 

In the critical, non-dissipative model, all mixed 2-point correlators of boundary height 
variables decay like r~ 4 . 



6.3. Other lattices 

The sandpile can be defined and studied on any type of lattice, but the most interesting 
examples are arguably the two-dimensional regular lattices. Among these, the triangular or 
the honeycomb lattices are natural to test the universality properties of the model. Little 
has been done so far. 

Numerical simulations have been carried out on the square, honeycomb and triangular 
lattices to determine the spectrum of exponents of toppling waves, which showed that the 
exponents are identical for the three types of lattices [78]. 

Correlations of heights 1 have been computed exactly on the honeycomb lattice in [79] . 
The results confirm that the height 1 scales in exactly the same way as on the square lattice, 
in the bulk and on an open or closed boundary. Only the lattice normalizations differ. 

The techniques developped recently in [35] to compute LERW passage probabilities offer 
a very interesting perspective since they can also be used on the triangular and honeycomb 
lattices. The results contained in [35] already yield the height distribution at one site. On 
the honeycomb lattice for instance, they read (each site in the bulk has three neighbours so 
that the heights take three values) 

Pi = — , P 2 = — , P 3 = -• (6.5) 
12' 24' 8 K ' 

The calculation of 2-point correlations involving a height 1 and a higher height should 

be possible using these techniques and would presumably confirm the type of scaling found 

on the square lattice, namely that height 1 is not logarithmic while the higher heights are. 

These computations remain to be done. 



6.4- Conformal invariance from SLE 

In recent years, the search of conformally invariant properties has been generalized to 
extended geometrical objects, like certain clusters and interfaces. In this respect the most 
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successful approach is the one based on the Stochastic Loewner Equation (SLE), which has 
led to an enormous body of results, see for instance the reviews [80, 81]. 

These ideas have been applied to the ASM on the square lattice in [82] by looking at the 
avalanche clusters. More precisely, after the addition of a grain, an avalanche may take place, 
during which a certain number of sites topple, usually several times. The set of sites which 
topple at least once form the avalanche cluster. By generating these clusters numerically, the 
geometric features of the avalanche cluster boundaries can be studied. 

It has been argued [82] that the avalanche frontiers are random curves with a fractal 
dimension numerically very close to 5/4, the same value as the loop erased random walks. In 
addition the statistics of these curves appears to be consistent with a description by an SLE 
process for the value of the parameter k — 2, corresponding to a value c = — 2 of the central 
charge. The traces of SLE K=2 are known to be statistically identical to the LERW [83]. 

Similar simulations on the honeycomb lattice [79] lead to the same conclusion, confirming 
its universality. 

It is important to note that the clusters examined here belong to the dynamical aspects 
of the ASM, since they appear during the relaxation process. Consequently they are not 
the kind of random variables whose statistics are in principle determined by the stationary 
measure. 

7. Conclusion and perspectives 

In this review, we have collected the most important results obtained in the stationary regime 
of the Abelian sandpile model, in connection with the issue of conformal invariance. The 
main results concern: the height variables, in the bulk and on boundaries, different types of 
boundary conditions and the fields effecting a change of boundary conditions, the crucial role 
of dissipation. A large number of mixed correlation functions have been calculated on the 
lattice, either analytically or numerically, and perfectly match the conformal predictions. 

The fairly safe conclusion we may draw from these data is that the ASM in its stationary 
regime is a conformal invariant system, and that it provides a lattice realization of a 
logarithmic conformal theory. We believe that it is one of the few lattice models where 
the logarithmic conformal invariance, with its strange and unusual features, can be best 
understood, and is a rare example where correlation functions can be explicitly computed. 

Our current understanding of the model from the conformal point of view remains 
nonetheless very lacunar. The main question of course is to improve our knowledge of the 
spectrum of representations present in the conformal theory, both bulk and boundary, and to 
continue to investigate their precise nature and their fusions. The task can be tough; we have 
seen no sign so far of an extended symmetry, meaning that the relevant conformal theory 
could very well be irrational. Even though certain observables are correctly accounted for 
within the free symplectic fermion theory, the conformal nature of the bulk height variables 
suggests that the latter is in fact not the correct theory to describe the physics of the ASM. 
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